Observable Signatures of EMRI Black Hole Binaries 
Embedded in Thin Accretion Disks 



O 

< 
(N 



< 

o 

6 



> 

(N 
(N 
m 
(N 

^' 
O 



Bence Kocsis/ Nicolas Yunes,^' ^ and Abraham Locb^ 

^Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA. 
^ Dept. of Physics and MIT Kavli Institute, Massachusetts Institute of Technology, 
77 Massachusetts Avenue, Cambridge, MA 02139, USA. 
(Dated: January 19, 2013) 

We examine the electromagnetic (EM) and gravitational wave (GW) signatures of stellar-mass 
compact objects (COs) spiraling into a supermassive black hole (extreme mass-ratio inspirals or 
EMRIs), embedded in a thin, radiation-pressure dominated, accretion disk. At large separations, 
the tidal effect of the secondary CO clears a gap. We show that the gap refills during the late GW- 
driven phase of the inspiral, leading to a sudden EM brightening of the source. The accretion disk 
leaves an imprint on the GW through its angular momentum exchange with the binary, the mass 
increase of the binary members due to accretion, and its gravity. We compute the disk-modified GWs 
both in an analytical Newtonian approximation and in a numerical effective-one-body approach. We 
find that disk-induced migration provides the dominant perturbation to the inspiral, with weaker 
effects from the mass accretion onto the CO and hydrodynamic drag. Depending on whether a gap is 
present, the perturbation of the GW phase is between 10 and 1000 radians per year, detectable with 
the future Laser Interferometer Space Antenna (LISA) at high significance. The Fourier transform 
of the disk-modified GW in the stationary phase approximation is sensitive to disk parameters 
with a frequency trend different from post-Newtonian vacuum corrections. Our results suggest that 
observations of EMRIs may place new sensitive constraints on the physics of accretion disks. 

PACS numbers: 04.30. Tv,98. 62. Mw,04.30.-w,95.30.Sf 
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I. INTRODUCTION 

The full exploitation of gravitational wave (GW) sig- 
nals will hinge on controlling all systematics associated 
with their astrophysical sources. One can classify such 
systematics into three major groups: instrumental, the- 
oretical, and astrophysical. Instrumental systematics 
are associated with possible issues related to the detec- 
tor. For example, the future Laser Interferometer Space 
Antenna (LISA) might suffer from instrumental 

glitches [5|. Such glitches, and other potential instru- 
mental issues, might lead to a foreground of noise arti- 
facts that might have to be either removed, or dealt with 
via data analysis techniques. 

Theoretical systematics are due to incomplete model- 
ing of waveform templates @. The extraction of GWs 
from noisy data requires the construction of these opti- 
mized filters, which represent our best guess model for 
the GWs generated by the source. Since approximation 
schemes (either analytical or numerical) are employed to 
solve Einstein's equations of General Relativity (GR) for 
the source, the templates used for data analysis are not 
exact solutions and can introduce errors in parameter es- 
timation @. 

Astrophysical systematics arise from modifications to 
the waveforms caused by the environment. For ex- 
ample, when modeling GWs from black hole (BH) or 
neutron star (NS) binary coalescences, one usually as- 
sumes the binary is isolated from external perturbers and 
ambient electromagnetic (EM) or matter fields. How- 
ever, the unresolved GW foreground of Galactic and ex- 
tragalactic white dwarfs (WDs), and possibly extreme 



mass ratio inspirals, introduces additional astrophysi- 
cal noise for LISA sources (see Ref. and references 
therein). Furthermore, an additional nearby supermas- 
sive BH (SMBH) in the vicinity of a merging binary can 
lead to detectable Doppler-shifts in the GW signal Q. 

Astrophysical systematics are expected to be negligi- 
ble for the final stages of the inspiral and merger of two 
SMBHs, because the SMBH's inertia greatly exceeds that 
of the environment. However, this is not the case for ex- 
treme mass-ratio inspirals (EMRIs), where a small com- 
pact object (CO) spirals into a SMBH In this case, 
the GW inspiral rate and signal amplitude are decreased 
by a factor of the mass-ratio, making these systems more 
sensitive to astrophysical perturbations as well as the- 
oretical uncertainties. EMRIs produce millions of GW 
cycles in the LISA frequency band with signal-to-noise 
ratios (SNRs) around 20 with a GW phasing accuracy 
better than 1 radian and logarithmic mass measurement 
accuracy of 10~'^-10~^ for a typical source at 1 Gpc ob- 
served for a year (loj . 

In this paper we consider the most important effects 
that a radiatively-efficient, thin accretion disk might have 
on an EMRI that is embedded in it: 

(i) SMBH mass increase due to accretion; 

(ii) CO mass increase due to accretion; 

(iii) modification of the gravitational potential due to 
the disk's self-gravity (e.g. changing the angular 
velocity of the orbit 0jS &- function of radius, and 
inducing additional apsidal and nodal precession); 

(iv) modification in the energy and angular momen- 
tum dissipation rate (e.g. hydrodynamic drag from 
winds, torques from spiral arms, and resonant in- 
teractions analogous to planetary migrations). 

We examine the conditions necessary for the tidal gravity 
of the CO to open a gap in radiation-pressure supported 
accretion disks, and its implications on the EM and GW 
signals. In particular, we study whether LISA will have 
sufficient sensitivity to resolve the presence and structure 
of the accretion disk. 



A. Relevance of Accretion Disks to EMRIs 

EMRIs are expected to form in dense galactic nuclei 
of stars, WDs, NSs, and BHs in orbit around a central 
SMBH. These dense nuclei arc sometimes called galactic 
cusps (see e.g. Q), because of their sharply peaked den- 
sity profile at zero radius. Some of these galactic nuclei 
are coincidentally active, meaning that gas is currently 
accreting onto the central SMBH and produces bright 
EM radiation. Accretion disk effects on EMRIs are most 
prominent in active galactic nuclei (AGN) where gas is 
actively feeding the central SMBH. 

Plausible arguments have been put forth both against 
and in favor of the common existence of EMRIs in AGN 
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disks. EMRIs can only be detected at relatively low red- 
shift (2 < 0.5) [ll|, but only a small fraction of galaxies 
within z < 2 are active and host a massive gaseous disk. 
AGN activity may be triggered by the inflow of gas dur- 
ing major galaxy mergers [l^l- However, the SMBHs in 
the centers of merging galaxies form a binary, which may 
deplete the central cusp of stars : thereby reducing the 
probability of EMRI events. On the other hand, AGN 
activity may be fueled by the tidal disruption of stars in 
dense central cusps, which have large EMRI rates (l^ . 
Stars may be captured or may form in accretion disks 
by fragmentation and/or coagulation of density enhance- 
ments The remnants of these stars would be 
pushed inwards by the disk and could provide a reservoir 
of EMRI events in AGNs. 

Astrophysical evidence already exists for tightly bound 
BH binaries with accretion disks. OJ287 is believed to 
be an SMBH-SMBH binary, with masses of IO^Mq and 
10^'^ Mq; respectively, orbiting in an inclined accre- 
tion disk of mass ~ IO^Mq [l^. For this system, opti- 
cal flashes are observed periodically and interpreted as 
crossings of the accretion disk by the smaller object. For 
EMRIs with stellar mass COs, similar EM flares will be 
much harder to detect. 

Estimates on the expected EMRI rates are very uncer- 
tain, around a few tens to a few hundreds per year (Tlj . 
LISA is expected to be sensitive to EMRIs up to redshifts 
close to unity Q , although most events should be at red- 
shift much smaller than unity. A few percent among these 
might be in AGN environments, where accretion disk ef- 
fects are non-negligible. These sources, if observed, may 
become the most interesting EMRI sources for studying 
astrophysics with LISA. 

We examine whether an EMRI, if present in the ac- 
cretion disk, is capable of regulating the accretion of the 
SMBH. If the secondary is sufficiently massive, its tidal 
gravity would expel gas from the inner regions, greatly re- 
ducing the amount of gas that would fall into the SMBH, 
and thus decreasing the disk's EM luminosity. However, 
as the inspiral proceeds, and the relative importance of 
tidal field changes compared to the local radiation pres- 
sure in the disk, the disk might refill, reigniting bright 
AGN activity. 

The presence of an accretion disk around EMRIs also 
leads to interesting possibilities for future GW detections. 
A detection of the imprint of an accretion disk on the GW 
signal could inform us about accretion disk physics. But 
by the same token, the presence of the disk complicates 
the modeling of the GW signal, as the rather uncertain 
accretion disk physics introduces additional theoretical 
errors and potentially makes tests of GR with EMRIs 
more difficult. 

If GWs could inform us about accretion disk physics, 
they could then provide candidates for EM counterpart 
searches. Indeed, accretion disks are very sensitive to 
the accretion rate parameter. A LISA measurement of 
the accretion rate would imply a plausible range of AGN 
luminosity. Then, by looking at the LISA source location 



box (approximately 1° angular, 10~^ distance measure- 
ment accuracy [lO|), EM instruments could search for 
AGNs with the predicted luminosity and redshift [l9j . 
Since bright AGNs are relatively sparsely distributed in 
the universe relative to the LISA error volume, a search 
could cut down the number of galaxy candidates to just 
one Alternatively, the EM counterpart could be 

identified in case it is strongly modulated in time by the 
EMRI. 

The identification of an EM counterpart would allow 
the use of EMRIs as standard sirens: measuring the 
distance using the EM redshift and the GW luminos- 
ity, which would allow to independently test cosmologi- 
cal models [2l| - [23j . Peculiar velocities and weak lensing 
are expected to be the main limitation for studying cos- 
mology with LISA, implying that a large number of low 
redshift sources such as the EMRIs considered here, are 
necessary for tightening existing cosmological constraints 



B. Previous Explorations 

Accretion disks are common astrophysical systems that 
have been studied in depth, but only recently has there 
been some effort to discuss their effect on GW sources. 
To our knowledge, the first study of accretion disk effects 
on GWs was by Giampieri He considered an equal- 
mass binary, where each component is gaining mass due 
to accretion, in turn leading to a modified radial inspiral 
rate. He then argued that a measurement of the so-called 
GW braking index, kaw = fcwfcvj/ fow^ where /qw is 
the GW frequency, could lead to information on the rate 
of accretion. However, if limited by the Eddington rat^B, 
the accretion timescale is ^ 10* times larger than the 
observation time and the effects on GWs is insignificant 
(see Sec. IIVAI below). 

Almost simultaneously with Giampieri, 
Chakrabarti |25|. [26| and then Giampieri, Gcrardi 
27[ considered an EMRI embedded 
disk, where the CO accretes in the 
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Bondi-Hoyle-Lyttleton (BHL) approximatioifl [281-133. 
Assuming the gas velocity is different from the CO's 
velocity due to pressure gradient effects, they found 
that the CO experiences a head-wind (accelerating the 
coalescence) or a tail-wind (delaying the coalescence), de- 
pending on whether the disk is rotating at sub-Keplerian 
or super-Keplerian velocities. For very high accretion 
rates, the flow is trans-sonic and super-Keplerian, and 



^ The Eddington limit is defined as the mass accretion rate at 
which the inward gravitational force is balanced by the radiation- 
pressure force produced by the in-falling matter in spherical sym- 
metry; see Sec. HVAl for details. 

^ BHL accretion results when the accreting object is completely 
embedded in a gaseous medium and accretes isotropically; see 
Sec. IIV Bl for further details. 



4 



the tailwind could supersede the angular momentum 
loss via GW emission, leading to a stalled orbit or even 
an out-spiral. 

Later on, Narayan [s^l examined the effeet of 
radiatively-ineffieient accretion flows on EMRIs. In qui- 
escent nuclei, the accretion rate is often much lower than 
the limiting Eddington value and the accretion flow is 
dominated by advection of the thermal energy. Accre- 
tion disks of this type are commonly referred to as be- 
ing "radiatively inefheient," since the thermal energy of 
the gas is advected inwards rather than radiated away 
as in thin accretion disks. Narayan estimated the impor- 
tance of hydrodynamic drag by computing the ratio of 
the timescale on which the EMRI loses angular momen- 
tum due to hydrodynamic torques and GW emission. For 
all reasonable sets of ADAF disk parameters, he found 
that the GW phase is changed by order 10~^ radians, 
well below the measurement accuracyof LISA. 

Shortly after, Subr and Karas H, investigated the 
evolution of the eccentricity and inclination of a CO when 
crossiiig a thin radiatively-efhcient Shakura-Sunyaev a- 
dislJl [Sg. Assuming that the disk is not perturbed by 
the CO significantly, they found that the orbit circular- 
izes and aligns with t he p lane of the disk in its outskirts. 
In a follow-up work [33, they extended their study to 
the inner parts of the Shakura-Sunyaev disks, where the 
radiation pressure significantly modifies the density pro- 
file, and examined hydrodynamic drag during disk cross- 
ing, and angular momentum exchange with spiral density 
waves analogous to Type I and II planetary migration. 
They have provided formulas for the relative timescale of 
the effects, and concluded that GW emission drives the 
evolution interior to ^ 100GA//c^. 

Another important study in this field was pre- 
sented by Levin (T^]. Motivated by simulations of 
Gammie [s^l, Levin constructed models of thin self- 
gravitating, radiatively-efhcient disks, including optically 
thick (i.e. photons scatter several times before escap- 
ing the disk) and thin regimes. He re-derived order- 
of-magnitude estimates of the non-relativistic hydrody- 
namic drag and planetary migration timescales, similar 
to that of Karas and Subr , and also included the ef- 
fect of azimuthal winds. Based on these estimates. Levin 
argued that, although GW emission drives the evolution 
inside ^ 100GA//c^, such disk effects may be important 
for LISA EMRIs. 

More recently, Barausse and Rezzolla [s^ examined 
the effects of relativistic hydrodynamic drag on EMRIs 
embedded in a thick torus. They found that the hydro- 
dynamic drag drives the EMRI toward alignment over- 
coming the GW radiation reaction, which by itself would 
drive the orbital plane toward anti-alignment with the 
MBHs spin [s^. However, stationary thick massive tori 
with constant specific angular momentum are unstable 



^ See Sec. Illll for further details on Shakura-Sunyaev Q-disk mod- 
els. 



to global non-axisymmetric modes that grow on a dy- 
namical timescale (40l - |4^ . making their conclusions on 
EMRIs rather uncertain. 

The EM emission of accretion disks around sub-parsec 
scale BH binaries have been considered by many au- 
thors. However, previous explorations focused on com- 
parable mass binaries. In this case, the gravitational ef- 
fects would clear a gap around the binary, significantly 
reducing accretion. As long as the binary separation is 
large enough that the gas can follow the GW inspiral 
rate of the binary, tidal stresses would act to increase 
the EM luminosity of the disk Eventually, however, 
the gas is left behind at radii exceeding lOOGM/c^ for 
equal mass systems and the gap freezes relative to the 
rapid GW-induced inspiral. Consequently, bright EM 
emission would be expected only several years after the 
merger, when the gas has had time to diffuse inwards and 
accrete onto the remnant SMBH However, periodic 
inflow across the gap may generate EM variability prior 
to merger (45l - l48| or the secondary can shepherd any gas 
remaining interior to its orbit into the SMBH after gap 
freez-out H^. The conditions for gap opening have not 
been examined in radiation-pressure dominated accretion 
disks, which is most relevant at separations smaller than 
^ 1000Gil//c^. Here, we investigate whether this can 
lead to gap refilling prior to merger, which could have 
important consequences for EM signatures of EMRIs. 

None of the previous studies examined in detail 
whether LISA has the sufficient accuracy to resolve the 
imprint of accretion disk effects on EMRI GWs and 
whether any of the accretion disk parameters may be 
recovered; this is the main topic of this paper and a 
companion paper [s^] (hereafter Paper I). In Paper I, 
we examined the detectability of the angular momentum 
exchange of EMRIs with an ambient accretion disk using 
a large class of torque models parameterized by two free 
parameters. In particular, we pointed out that plane- 
tary migration models are examples that generate a very 
significant GW phase shift for LISA. 

In this paper, we begin by reviewing the astrophysi- 
cal models of the most important accretion-disk effects 
on EMRIs in AGNs, many of which were not included 
in Paper I. We focus on standard radiatively-efficient 
thin disks, where the viscosity is proportional to the to- 
tal pressure (Shakura-Sunyaev a-disks [s^) or propor- 
tional to gas pressure only (which we refer to as /3- 
disks) [5l|. These disks constitute the standard model 
of luminous AGN accretion disks (see Refs. [3l|, [52j 
and Sec. IIIII below). For both of these models, we in- 
clude a more detailed analysis of the effects considered 
previously (SMBH mass accretion, CO mass accretion, 
hydrodynamic drag, torques from spiral density waves 
and resonant interactions), investigating the detailed cir- 
cumstances under which such effects are possible or sup- 
pressed, and derive the effects of axisymmetric disk self- 
gravity on GWs, which was not considered before. We 
provide a detailed study of the consequences on the GW 
observable, using a few GW data analysis tools. 
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The analysis presented here is by no means exhaus- 
tive. For simphcity, we restrict attention to EMRIs on 
non-inchned, quasi-circular orbits. This restriction does 
not allow us to model binaries in which eccentricity is ex- 
cited in the presence of an accretion disk [4^, [ZJ [ssl - lssj . 
Other effects that we do not discuss here include the fol- 
lowing: GWs generated by accretion flows through the 
excitation of BH quasinormal modes [56l - (6l| : GW en- 
ergy flux dissipation by an ambient viscous disk, driving 
transverse and longitudinal density waves [G^I that could 
heat the dis k sig nificantly and result in an observable in- 
frared flare |63|; EM radiation generated by GWs in a 
strongly magnetized plasma, boosting the frequency of 
photons [6^l65l|. driving magnetosonic waves [Gq . focus- 
ing of EM radiation |67ll , generating photons in a static 
magnetic field [gsI . |69| and photons back-converting to 
GWs in a magnetized plasma (70l - [7^ . We also do not 
consider the direct scattering of GWs by the gravity of 
the gas, nor the GW radiation of spiral arms in the disk. 
Clearly, these effects are interesting and should be stud- 
ied in more detail, but they go beyond the scope of this 
paper. 



C. Executive Summary of Results 

This sub-section of the introduction is an executive 
summary of our main results, intended for non-specialists 
or for readers who might not be interested in all the re- 
lated technical details of the paper, considering its length. 

We derive the necessary conditions for the tidal ef- 
fect of the CO to open a gap in radiation-pressure dom- 
inated a and /3-disks. We find that EMRIs can open 
gaps at large radii in both a and /3-disks. Depend- 
ing on the EMRI masses and accretion disk parame- 
ters, the gap typically closes during the inspiral due to 
strong radiation-pressure gradients. Gap refilling occurs 
at orbital separations outside (inside) the LISA frequency 
band at orbital radii (in geometric units, sec Scc. lID[) r > 
300 (24 M.) for a (/3) disks, under typical EMRI parame- 
ters [SMBH and CO masses of (M,,m^) = (10^ IO)Mo]. 
Complete gap refilling occurs within 9 months in /3-disks 
for these EMRI masses, well before the inspiral termi- 
nates in coalescence. This implies that bright AGN ac- 
tivity is coincident with LISA EMRIs in both cases. 

We calculate the perturbations of the GW waveforms 
for quasi-circular, non-inclined EMRI orbits, due to the 
effects of mass accretion onto the SMBH and the CO, 
the hydrodynamic drag caused by an azimuthal and ra- 
dial wind, and the axisymmctric and non-axisymmetric 
gravitational effects of the disk. Some of these effects 
lead to a strong imprint on the GW phase, while others 
do not. 

We find that the gravitational torque from spiral den- 
sity waves in the accretion disk (also known as migration 
torque in planetary dynamics) provides the most signifi- 
cant deviation from vacuum EMRI dynamics. Migration 
leads to the strongest GW imprints when the CO clears 



a gap, leading to gas accumulation near the outer edge of 
the gap, like in a hydroelectric dam. The CO's mass ac- 
cretion also leaves a significant imprint on the GW signal, 
if described by spherical BHL accretion [l3l. However, 
we find that the CO accretion rate is significantly re- 
duced by many processes, predominantly by limited gas 
supply and radiation pressure. Accounting for these lim- 
itations, the GW signature due to the CO's mass accre- 
tion is typically much smaller than that of migration. If 
tidal torques from the CO open a gap, the mass accretion 
onto the CO may be greatly reduced. All other effects 
are typically less significant. 

We calculate the effect of all of the above mentioned 
processes on the GW observables, namely the waveform 
amplitude and phase, both with a leading-order New- 
tonian waveform model as well as with a relativistic 
effective-one-body (EOB) model f73l - [75| . Figure [U sum- 
marizes the imprint of such processes on the Newtonian 
phase for the dominant GW mode as a function of the 
final orbital radius for a one-year observation. Different 
curve colors correspond to the phase difference between 
vacuum GW phases and those that include various disk 
effects, while solid and dotted lines correspond to a and 
/3-disks respectively. We provide simple asymptotic an- 
alytical formulas describing the phase shift perturbation 
for arbitrary accretion disk, EMRI, and observation time 
parameters in Eqs. (|TTCHTTB)) and Table M in Sec. IVlIIl 
Figure [T] shows that the GW phase is modified signifi- 
cantly for EMRIs relevant to LISA (those with masses 
(A/., TO*) = (I0^ 10) and final orbital radii r/ < 50M.) 
at 1 Gpc and an observing time of one year. The upper 
and lower thick magenta lines represent a rough measure 
of LISA sensitivity to the phase shift for a source at 1 Gpc 
and lOMpc, respectively. 

The various curves in Figure [T] exhibit interesting fea- 
tures at different radii, which correspond to different as- 
trophysical mechanisms that come into play. Most no- 
tably, the big decrease for the blue migration curves at 
r < 24 M, correspond to a transition from Type-II to 
Type-I migration as the gap refills for /3-disks. Coinci- 
dcntally, roughly interior to that radius, BHL accretion 
and hydrodynamic drag from azimuthal winds are acti- 
vated. The wiggles and the cutoff in the BHL accretion 
induced phase shift (black curves) correspond to a variety 
of effects. For a-disks (solid black curve), the gas den- 
sity and sound speed determine the BHL accretion rate 
at small separations. At larger radii, the Bondi accretion 
radius becomes larger than the disk thickness, reducing 
the accretion rate. Even farther out, differential rota- 
tion of the disk reduces the BHL rate and the decrease 
in the average background radial gas velocity makes the 
amount of gas supply an important limitation. In the in- 
nermost region of the disk, photon diffusion is slow and 
the radiation is trapped within the BHL fiow leading to 
super-Eddington accretion rates. However, at r > 35A'/, 
this is no longer true for an a-disk and the flow becomes 
Eddington limited, greatly reducing BHL accretion ef- 
fects. 
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FIG. 1. The GW phase shift as a function of final radius in units of A/, induced by different accretion disk effects relative to 
vacuum waveforms (see Sec. II DI for the conventions used here). Solid (dotted) curves correspond to a. (/3) disks, with different 
colors indicating different disk effects: black corresponds to Bondi-Hoyle-Lyttleton (BHL) accretion, green to azimuthal wind 
and blue to migration. The thin, solid magenta line is the total accumulated GW phase in vacuum. The thick, solid (dashed) 
magenta line corresponds to a measure of the accuracy to which LISA can measure the GW phase for a source at f Gpc 
(10 Mpc). Observe that certain disk effects, like migration, can leave huge imprints on the GW observable, inside the LISA 
accuracy bucket. 



Relativistic waveform models yield similar results to 
those presented in Fig. [T] After aligning the wave- 
forms in time and phase (equivalent to a maximization 
of the SNR over the corresponding extrinsic parameters 
in white noise), we find changes in the GW phase after 
a typical one-year inspiral of up to 0(10'*) radians when 
modeUng migration in /3-disks. Migration effects for a- 
disks are much smaller, since these disks are less dense. 
As the gap is typically expected to close for EMRIs in the 
most sensitive LISA frequency band (r < 25A/,), Type-I 
migration is the most relevant process. Supply-limited, 
BHL accretion and wind effect lead to a dephasing of 
0(1) rads. Other effects are less significant: 0(10""^) ra- 
dians for SMBH mass accretion, and 0(10^*) radians for 
axisymmetric self-gravity effects. 

We then proceed with a more careful data analysis 
study on the distinguishability of accretion disk effects by 
computing a data analysis measure for two representative 
systems at 1 Gpc with component masses (10, 10^) Mq 
and (10, 1O^)M0, respectively. We calculate the SNR 
in the waveform difference between signals accounting 
for accretion disk perturbations and those that do not, 
marginalizing over an overall time and phase shift. We 
find that, for these systems, p{5h) > 10 after just 4 
months of evolution for /3-disk migration, while it takes 
one full year of integration to reach the same SNR for 
BHL accretion and wind effects. All other accretion disk 
effects are less significant. 

Finally, we examine possible degeneracies between ac- 
cretion disk effects and vacuum EMRI parameters. We 
analytically derive the Fourier transform of the wave- 
forms in the stationary phase approximation. We find 
that the disk-induced perturbation to the frequency- 



domain GW phase depends on the GW frequency to 
a high negative power relative to the Newtonian term, 
multiplied by a function of the initial binary masses, 
the a-disk parameter and the SMBH accretion rate m,. 
In contrast, the phase of the Fourier transform of vac- 
uum waveforms is a positive power of frequency relative 
to the Newtonian term, when including post-Newtonian 
(PN) corrections. The difference in the frequency scaling 
arises because the accretion disk effects grow with orbital 
separation (lower frequency), as opposed to PN correc- 
tions which grow with decreasing separation (higher fre- 
quency). This suggests that accretion disk effects are 
not strongly correlated with general relativistic vacuum 
terms in the frequency-domain GW phase. Whether this 
statement holds in a realistic data analysis implementa- 
tion requires a much more detailed analysis of the likeli- 
hood surface that is beyond the scope of this paper. 

Our results suggest that if a GW signal is detected 
from an EMRI in an accretion disk, then matched fil- 
tering with accretion disk templates could allow for the 
measurement of certain disk parameters to an interesting 
fractional accuracy (early estimates suggest 10% accu- 
racy for certain parameters [76|). The precise magnitude 
of the latter requires the detailed mapping of the likeli- 
hood surface with relativistic EMRI signals, full Fourier 
transforms, and improved disk modeling (including rela- 
tivistic effects and magnetic fields), which is beyond the 
scope of this paper. 

We caution, however, that the models considered here 
might not provide a fully realistic description of the an- 
gular momentum exchange between the binary and the 
accretion disk, leading to systematic theoretical uncer- 
tainties in interpreting GW measurements. Since EM 
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observations of the disk luminosity are also sensitive to 
m,, the combination of contemporaneous EM and GW 
observations might hold the key for constraining the ac- 
cretion disk physics most reliably. 



provides a more detailed analysis that includes higher- 
order rclativistic effects. 



A. Basics of EMRI Dynamics 



D. Organization and Conventions 

This paper is aimed at both the General Relativity 
and Astrophysics communities. We present a significant 
amount of background material to make the paper self- 
contained for both communities. Sec. HIl reviews the ba- 
sics of EMRIs as relevant to GW physics and a rough 
measure of the accuracy to which the waveforms need to 
be computed for LISA parameter estimation. Sec. IIIII 
presents the basic elements of the thin accretion disk 
models that are considered in this paper, and derives 
the necessary conditions for tidal effects to open a gap in 
the disk around the secondary. Sec. IIVI studies the effect 
of binary mass increase due to gas accretion on the GW 
signal. Sec.|V]focuscs on the effect of hydrodynamic drag 
on the GW signal, induced by the gas velocity relative 
to the GO (i.e. wind). Sec. |Vl] discusses the effects of the 
axisymmetric gravity of the disk. Sec. IVIII concentrates 
on gravitational angular momentum exchange with the 
disk (i.e. migration) and its effects on the GW signal. 
Sec. IVIIII compares and contrasts the effect of the dif- 
ferent accretion disk effects on the GW phase. Sec. IIXI 
describes the theoretical framework through which we 
compute effective-one-body waveforms in the presence of 
an accretion disk. Sec. |X] performs a simple data analysis 
study to infer the detectability of accretion disk effects. 
Finally, Sec. IXII concludes and suggests future work. 

Throughout the paper, we employ the following con- 
ventions. We use geometric units with G = c = 1 un- 
less otherwise noted. This implies that masses arc in 
units of length or time, where the mapping is simply 
Mq = 1.476 km = 4.92 ^s. The EMRI is assumed 
to be composed of a SMBH with mass M, and a GO 
with mass m^. The SMBH is assumed to be spinning 
with spin angular momentum S, = a,M,, aligned with 
the orbital one. We do not model here the spin of the 
GO. We measure quantities relative to their typical mag- 
nitudes and denote Aj, = A/{10^Mq). For example, 
A/.5 = M,/{10^Mq) and m^i = m*/(10MQ). The ra- 
dial orbital separation is always scaled in terms of the 
SMBH's mass, such that f = r/M,. The natural scale 
for the start and end of observation in the most sensi- 
tive part of the LISA band is 20M, and lOM,, so we use 
^20 = r/[20M,) and fio accordingly. We use r' to denote 
distance from the GO to a field point. 



II. REVIEW OF EMRI GWS 

We start by reviewing some basic facts about EMRI 
dynamics, focusing only on leading-order effects. Sec. IIXI 



Since the GO orbits very close to the SMBH, GR ef- 
fects cause the largest perturbations of Newtonian orbits. 
In this section we consider eccentric EMRI dynamics, al- 
though in most of what follows we restrict our attention 
to quasi-circular EMRIs. In the absence of spin, the bi- 
nary's energy is 



2 r 2r 

the Kcplcrian orbital frequency is 



(1) 



(2) 



where f = r/M, is the semi- major axis of the orbit in 
units of A'/, . Apsidal or pericenter precession for an orbit 
with eccentricity e is 



GR,ap 



= 3A/.-i(l-e2)-if-5/2, 



(3) 



while Lcnsc-Thirring precession of the node of an inclined 
orbit around a spinning SMBH is 



f^LT = 2a, M, 



1^-3 



(4) 



The GO inspiral produces GWs that remove binding 
energy and specific angular momentum from the system 
at the rate 

i?ow ^ ^ ^ = -6.4 X 10-^^ ^ ^,(5) 



5 M2 f5 

32 g2(e) 

' 5 M, f7/2 



m2 =5 
-""•S '10 



to leading order in m^, <C A/,, where rio = 7'/10 and 
where 



51 (e) 



1 + Z3e2 , 37 4 
(1 - e2)7/2 



32(e) 



8^ 



(l-e2)^ 



• (7) 



We have here introduced the notation Ab = A/(\Q^y[Q) 
for any quantity A, such that Af.s = Af,/(10^Mo) and 
m^i = m^/(10MQ). 

Such loss of energy and angular momentum leads to 
the decrease of the semi-major axis and eccentricity at a 
rate 



64 TO* __ 



. 304 

"^<^-" = T^Ai^'^^(^)" 

where we have defined 



-1.3xl0'« ^g,(e)f-3 



2 X lO^'-'Af, 



A^.5 



-12,,r-l"^*l 



(8) 



© 



Mi 



33 (e) 



(1 - e2)5/2 



10 



(9) 
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Equivalcntly, we can parameterize the orbit in terms of 
the change in the orbital frequency 



• 96 N— 11/2 



(10) 



For quasi-circular orbits, several formula simplify. For 
example, one can easily see that -Eqw = f^^Gw- Similarly, 
the inward inspiral velocity is simply given by Eq. ([5]) 



with gi(e) 
axis, Eq. (H 
orbits, 



1. The orbital evolution of the semi- major 
can be integrated [t^, [t^ for quasi-circular 




and where fg and ff are the initial and final separations 
in units of M, for an observation time T. Let us define 
the critical radius and observation time where f/ starts 
to deviate significantly from fo as 



^/,crit 
^crit 



5 

256 TO* 



24TOlfAfr//^TV^ 



0.031 yr^ 



M. 



5^:4 



'/,io 



(12) 
(13) 



The measured GW phase to leading order is twice the 
orbital phase for a quasi-circular orbit 0gw = 2(/)orb- The 
total accumulated phase for a quasi-circular inspiral is 
then 



= 2 



n{t)dt = 2 



16 m^, ^ 




(14) 



where r is the dimensionless observation time defined in 
Eq. ([Tl]) . and f/ is the final radius at the end of the ob- 
servation. Depending on whether the observation time is 
short or long compared to the inspiral timescale T/f^ <C 1 
or r/fi ^ 1, Eq. ([Ti|) becomes 



6gw'' « 4 X 10^ rads — ^ 



r -3/2 



1-6 



T- 



5/8 



C^«2xl0^ads ^ 



1-0.1 



,.5/4-5/2 ^ 
™5/8 ^5/8 



(16) 



for short observations or widely separated binaries (T <C 
Tci-it and ff ^ ^/,crit)j a-nd long observations or close-in 
orbits (T ^ Tcrit and r/ <C r/ crit), respectively!! 



4 We note that </>g{^^ is the standard PN expression for the phase 
evolution as a function of time [TqII . when the phase evolution 
culminates in merger. Individual EMRIs, however, may outlive 
LISA observations, which is why we choose to use Eq. I I14II along 
with the asymptotes Eq. I ll5tfl6ll . 



Equations (|15f[TB|) show that initially the GW phase 
accumulates quickly as ciT-\-C2T'^, but then saturates to 
a rate CsT^/*, where Ci are time- independent constants. 
Setting TO,^ = IOM0 and M, = IO^Mq, we find that 
(/)gw ~ 0(10^) rads in a one year observation. 



B. Measures of LISA Sensitivity to GWs 

GW detectors are most sensitive to the phase of the 
GW signal. Since EMRIs can accumulate millions of GW 
cycles in the detector's sensitivity band, they make for 
excellent probes of the astrophysical environment. 



1. Simple Mass and Time-Scale Measures 

A rough measure of the accuracy to which LISA can 
extract parameters can be derived by looking at the mass 
measurement accuracy. This quantity is of order [loj 



(5M, SiTLi^ 

Alt TTli, 



10 



-3 



(17) 



where p is the SNR (see below). The SNR can be roughly 
as large as a few tens for a 10 M© CO spiraling into a 
10^ Mq within f < lOAf at a distance of 1 Gpc. Thus, 
the relative mass estimation precision is at best around 
10"** — 10~^ [l3|- Of course, for less distant sources, the 
SNR can be larger, allowing a better determination of the 
masses. This mass accuracy, compared to the accretion 
mass or local disk mass, provides a rough measure of 
whether the disk may generate important perturbations 
for LISA. 

Another rough measure to decide whether certain ef- 
fects are important for LISA is the following. If a pertur- 
bation has an associated timescale T on which it changes 
the inspiral phase by a factor of order unity, then the 
magnitude of the phase correction corresponding to this 
process is roughly 



r 



-tot 

''gw 



(18) 



where is the total, accumulated GW phase in the ob- 
servation. We assume here that this S(j) is not a simple, 
constant phase shift, but a modification in the phase evo- 
lution, such that at the end of the observation, the tem- 
plate dephases from the signal by an amount 5 (/)gw • LISA 
can detect phase differences of order a few radians (see 
below). Thus an effect is important if T < 4.4 x 10^ yrs 
for a single year observation, see Eq. p^ . In practice, 
T = x/i may be used where x{t) represents any of the 
following physical quantities: the mass of the SMBH or 
CO, the induced angular momentum or energy dissipa- 
tion rate relative to the GW-driven dissipation rates, the 
inspiral rate relative to the GW inspiral rate, or the fre- 
quency shift relative to the Kcplerian frequency. 
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2. Dephasing Measure 



Degeneracies and Template Placement 



A more accurate criterion to decide whether a certain 
GW modification is detectable for LISA is: 



> 



10/p ifp>10, 
if p < 10 , 



(19) 



where is the square of the signal-to- noise ratio (SNR), 
defined as 



df 



SnU) 



(20) 



The quantity h is the Fourier transform of the measured 
GW strain amplitude and S'„(/) is the spectral noise den- 
sity curve of the detector. Accounting for white-dwarf 
confusion noise and not averaging over sky-angles, the 
noise curve is 



Sn{f) = min| 



NSA g4.5Ty 



^N' ^NSA 



(21) 

where (S'^^SA ^.gai^ ^-cx-gai^ ^/-j j^^g functions of fre- 
quency, that can be found, for example, in Q- We 
account for sky position and binary orientation aver- 
aged response functions and noise curves, by multiplying 
Eq. dUl) by a prefactor of 20/3. 

Equation is motivated by the following arguments. 
Any accretion disk effect is measurable only if the EMRI 
is detected in the first place. We have here conservatively 
chosen p > 10 as the threshold for detection. Once the 
EMRI is detected, the accuracy to which a phase differ- 
ence can be measured is roughly 10/p, where p is the 
total SNR. 

Let us now relate this phase shift measure to the SNR 
of the waveform difference between signals that include 
and those that neglect disk effects: 5h= hi — h2- If these 
waveforms differ only in phase by an amount S(, 



6h = /ii (1 - e''^'^'^") 

The SNR of the difference is then 
df 



then 



(22) 



^n(/) 

d! 

Snif) 

df 

Snif) 



\Sh\\ 

\hi\^ [1 - cos((50Gw)] , 
Ihl^Scl^l^, (23) 



where in the last line we have assumed that Sc/jq^ <C 
1 rad. 

The perturbation of the waveform is significant if 
p{Sh) > 10, which is similar to the accuracy require- 
ments constructed in [sol - fs^ . If the instantaneous SNR, 
\hi\'^S~^{f) does not vary greatly while the phase differ- 
ence accumulates, then p{Sh) ^ {S4>Gw)p{hi) for (5</)gw <C 
1 rad, which leads to the simple phase shift criterion, 
Eq. HI]), above. 



One can generalize the above measures by allowing for 
both amplitude and phase modifications. The SNR of 
the waveform difference then becomes 



'{Sh) = min 
A2 



df 



Snif) 



h,if)-h2if;X2) 



, (24) 



where hi and /i2 are the Fourier transforms of two wave- 
forms (the "signal" and "template"), normalized such 
that pQii) = p(/i2) = 1- The template may depend on 
free parameters A2 that may be different from the true 
astrophysical values, where the minimum difference cor- 
responds to the best fit. Expanding Eq. ((24| . we find 
p^iSh) = 2MM, where MM = 1 - C is the mismatch 
and 



0{hi, ^,2) = max 
A2 



45R 



^n(/) 



hi{f)Kif;X2) 



(25) 



is the overlap. 

For a simple estimate, we minimize p'^{dh) over cer- 
tain non-physical extrinsic parameters only, i.e. an overall 
phase and time shift, as opposed to intrinsic parameters, 
such as the binary's masses or spins, or other extrinsic 
parameters, such as the distance to the source, the polar- 
ization angles or the sky position. When the SNR of the 
difference is computed in this way, estimates of measur- 
ability are optimistic as they do not account for possible 
degeneracies between all parameters. For example, a sig- 
nal hi{f, Ai) with true parameters Ai may be mimicked 
by a waveform model /i2(/, A2) with different intrinsic pa- 
rameters A2 7^ Ai , even if the extrinsic parameters are the 
same. In Sec. IX Cl we show that the incorrect determina- 
tion of EMRI parameters in a fiducial model (e.g. with 
no accretion disk) cannot mimic the waveform of a differ- 
ent model (e.g. including the effects of an accretion disk), 
because of the particular spectral features introduced by 
accretion disks in the Fourier transform of the response 
function. 



III. REVIEW OF ACCRETION DISK MODELS 

We proceed with an overview of the accretion disk 
models under consideration. As this paper aims to bridge 
between the accretion disk astrophysics community and 
GW physics community, we have chosen to provide a 
complete description of background material. 

A. Basics of Accretion-Disk Models 

Despite a long history of observational, theoretical, and 
numerical investigations, accretion disks remain one of 
the most exciting unsolved problems in astrophysics. The 
complexity is related to the modeling of magnetohydro- 
dynamic (MHD) flows, turbulence, radiative transport. 
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and plasma physics. Here we provide a brief overview 
of the formulas used to model accretion disks; for more 



details see textbooks by Shapiro & Teukolsky [31| and 
Frank et. al. js^ . 

We restrict our attention to geometrically thin, radia- 
tively efficient, stationary accretion disks, responsible for 
the observed bright emission around AGNs. In this case, 
as gas orbits around the central object, it radiates ther- 
mal energy away much faster than the timescale over 
which the gas particles drift inward, and the disk main- 
tains a thin configuration. Radiatively efficient disks are 
the most massive among accretion disks, as the mass ac- 
cretion rate is largest and the inward drift velocity is 
smallest. For lower radiative efficiencies around quiescent 
SMBHs, accretion is described by other models, such as 
advection-domination, which we ignore here. 

Radiatively-efficient, stationary accretion disks can be 
described by the Shakura-Sunyaev a-disk model [sl], [11] . 
The component of the viscous stress tensor corre- 
sponds to an effective viscositj[3 t^- = pvViVj and it 
is responsible for the slow inflow of gas. In the a- 
disk model, the viscous stress is assumed to be propor- 
tional to the total pressure Ptot in the disk at each ra- 
dius: = — (3/2)a Ptof The total pressure includes 
both thermal gas pressure and radiation pressure, and 
the dimensionless constant of proportionality a is a free 
model parameter. These disks are viscously, thermally, 
and convectively unstable to linear perturbations [83l - l86l | . 
In the alternative model [5l|, hereafter denoted /3-disks, 
viscous stress is proportional to the gas pressure only, 
tr(i> = — (3/2)q; Pgas, and such models are stable[f[ 

The nature of viscosity, however, is not sufficiently 
well-understood to predict which of these prescriptions 
is closer to reality. Recent MHD numerical simulations 
of accretion disks indicate that stresses correlate with to- 
tal pressure as in the a-disk model [s^] and are thermally 
stable [11] , though they might be viscously unstable [89j |. 
In cases where the diffusion scale is larger than the wave- 
length of the magneto-rotational instability, other simu- 
lations arc consistent with the /3-disk model [9^ [9l| . The 
instability in a-disks implies spectral variations which are 
not observed in many systems, while the /3-disk model 
provides a better match to spectral constraints [o^]. In 
this paper, we remain agnostic about the disk model and 
carry out calculations for both of them. 

Shakura-Sunyaev a and /3-disks for a BH of fixed M, 
mass are described by two free parameters: the accretion 
rate M, and the a parameter in the viscosity prescrip- 
tion. AGN observations show that the accretion rate rel- 
ative to the Eddington rat6[3, m, = M,/M,Edd, is typ- 



^ Throughout this section p = p{r) denotes gas density (i.e. not 
the SNR of GWs), v is the kinematic viscosity coefficient in units 
of cm^/s, and v is the gas velocity. 

^ Another popular model that is stable assumes t. 



gasPtot . 

^ We define the Eddington rate A/.e^j more precisely in Eq. 1541 1 
below. 



ically around 0.1-1 with a statistical increase towards 
higher luminosities [H, [13] . Theoretical limits based on 
simulations of MHD turbulence around BHs are inconclu- 
sive, but are consistent with a in the range 0.001-1 [osj . 
Its estimated value in protoplanetary accretion disks is 
lower, a < 10"'^ [9^. Observations of outbursts in bi- 
naries with an accreting WD, NS, or stellar BH imply 
a = 0.2-0.4 [H i98|. The value of a in AGN accretion 
disks is uncertain, but might be expected to be similar 
(for a review, see Ref. [l^). In the following, we assume 



m.i 



oT 



1, 



(1 



= 1, a-i 

0.1 



01 = ^' 



but retain m,i and ai to be able to describe different 
values. Here e is the radiation efficiency (see Eq. (jSTj) 
below). E 

The surface density, E(f), and the scale height of the 
disk from the mid-plane, H{f), can be calculated as [l6| . 



S(r-) 



H{f) 



37r3/5/2Kl/5 

27rc(l - /3)' 



kna 



4/5 



M.'/'l72/^/34(i-^)/5, 
(27) 

(28) 



where 6 = for a-disks and 1 for /3-disks, and the ra- 
dial dependence is implicit in the orbital velocity O = 
^-1^-3/2 g^j^^ Reie, 13(f) = PgasMot, where ptot = 
Pgas +Prad, Pgas is the thermal gas pressure, Prad is the 
radiation pressure, and /3 satisfies 



(l/2) + (l/10)(h-l) 



23/5^4/5^ / 



1-/3 



2/5 



(29) 



with /io = 2/{?,Xh + 1) = 0.615 the mean molecu- 
lar weight, fee the Boltzmann constant, ctsb the Stefan- 
Boltzmann constant, and /t = 3/4 a constant related 
to the assumption of optical depth (see [1^). The quan- 
tity K = /iedT/mp is the electron scattering opacity in 
the medium, with the number density of electrons 
relative to the total density, nip the proton mass, and 
ctt the Thompson scattering cross-section. In practice 
/ie = (1 + Xh)/2 = 0.875 and n = 0.348 cm^g-i for a 
fully ionized gas of hydrogen and helium, where the mass 
fraction of hydrogen is Xh = 0.75. Equation (PH)) can be 
solved numerically for /3(f) at each radius and substituted 
in Eqs. (|27H28l) to obtain E(f) and H{f). The kinematic 
viscosity coefficient in the disk is 



37rl] 



(30) 



* Note that in all of our formulas m,i and ei always appear in 
the combination rn,\/e\. To simplify notation we suppress the 
ei scaling. 
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[log . 

Equations (|27ff^ give a self-consistent non- 
relativistic description of a geometrically thin, ra- 
diativcly efficient, stationary accretion disk, provided 
the following conditions are satisfied: the disk is opti- 
cally thick to electron scattering, opacity is dominated 
by electron scattering, the disk is hot enough to be 
fully ionized, the self-gravity of the gas is negligible 
relative to the gravity of the accreting object, and the 
modifications near the inner boundary condition are 
neglected (see below). These conditions arc all satisfied 
within 6 < f < 10^ for M, < 10^ Mq, m > 0.1, a > 0.1 
jogf . Equation shows that well within 



cnn • 16/21 2/21 ,,r2/21 

^rad = 600 m,/ a/ M,^ 



(31) 



/? <C 1 holds, so that radiation pressure dominates over 
thermal gas pressure. In this case, Eq. (|29p can be in- 
verted analytically 

I3{f) « 4.3 X IQ-^a^^^^m,^^^ M,^^^^fll^^° for /3-disks 

(32) 

and Eqs. 
S(f). 



simplify to 
538.5 gcm^^ Qj^^m^j^f^Q^ for a-disks, 



-2 -1 • -1-3/2 

,6 „ „,^-2 ^,-4/5^3/5 ^^1/5 --3/5 

for /3-disks, 



1.262 X 10^gcni'^a^^'m,{°M,l r^Q 



^ ^™ i^'"io^ for a-disks. 



5.907 X 10-21 Mq -a^-m.^r. 

-1 „ -4/5 ■ 3/5 , .1/5--3/5 



1.384 X 10-^' Mo- 
tor /3-disks, 
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H{f) « 1.5 m.i M, = 1.5 X IO^Mqto.iA//. 



5 • 



(33) 
(34) 



The disk scale height is the same for the two models in 
the radiation pressure dominated regime, approximately 
constant in radius. Note that the thin disk assumption 
2H ^ r breaks down within f < 3m, i. Since radia- 
tively eflacient thick disks have no widely accepted an- 
alytical models to date, we extrapolate Eqs. p7l - [28l) to 
this regime, as well. 

Equation p3p shows that /3-disks are much more mas- 
sive within f ^ 1000. This is to be expected, as the 
effective viscosity is much smaller for /3-disks relative to 
a-disks by a factor of order /3. Since the orbital veloc- 
ity of gas is different for different radii, viscosity leads to 
energy dissipation, and a slow radial inflow. A smaller 
viscosity implies a smaller radial inflow velocity, which 
for a fixed accretion rate corresponds to a larger mass. 

The local disk mass near the binary in a logarithmic 
radius bin is defined as 



7.4 X 10^8 Mq 
for a-disks, 
1.7xlO-4Moar^/^m^fM:5^/^-[/^ 



(35) 



for /3-disks . 




1000 



FIG. 2. The velocity scales in the problem as a function 
of radius. From top to bottom (with corresponding equa- 
tions): Keplerian orbital velocity ([2]), differential rotation at 
the Hill's radius for fi- and a-disks (|37p . isothermal sound 
speed (|38p . azimuthal wind (|79p . radial inflow velocity for 
a~ and /?-disks (|37p . and the GW radiation reaction inspi- 
ral velocity (|5]). The scalings with accretion disk param- 
eters (ai,rh.i) are as labelled. Fiducial parameters used: 
ai = m.i = M,5 = m*i = 1. 



The gas density in the disk, the mean radial inflow ve- 
locity, and the isothermal sound speed are 



p(r) — 



(36) 



2.0 X 10-^^MQ-'^a^^ih,^M.-r 



1^^3/2 
10 



for a-disks 



2 „-4/5 . -2/5, .-4/5 -3/5 
1 "t.i JlJ. 'lo 



4.6x1O~23M0 -a, ' m.-, ' ' r 



wf"(f) 



for /3-disks 

3 z/ _ M, _ M,r 

2 r 2'KrYj 



(37) 



4„,.„;;-,2 ^ 5/2 Q,_disks 



-3.8 X 10 aim,;^r^(, 
-6.9 X lO-^ai/^m^/^Af^/Z^f^o^/^for /3-disks 



Cs(f) = y'^ = iJO w 0.047 m.ifj"(f/^ 



(38) 



where M, is the SMBH accretion rate, defined more pre- 
cisely in Sec. IIV Al 

Figure [2] compares the typical velocity scales in the 
problem for quasi-circular EMRIs. Typically, \vri,\ 

\Vr.p\ ^ \Vr,a\ ^ Cs for f > 10, whorC Vr.a and Vr.JS 

corresponds to vf'^^ for a and /3-disks. The figure also 
depicts other relevant velocity scales which we shall de- 
rive in Sees. |IV] and |V] below. Similarly, to get a feel 
of the typical disk mass scales. Figure [3] compares the 
MBH and EMRI mass, and the local disk masses for a 
and /3-disks, as well as the CO accretion rate per year 
derived in Sec. IIV Bl below. The local /3-disk mass and 



12 



10'" ^ 
10* - 
100 - 



M. 




§^ 


1 






ass 


0.01 






B 




- AMf 







lO"* 






lO"'' 










- Am^ 






10-8 




/ 

/ ^ 




10 



100 



1000 



FIG. 3. The mass scales in the problem. Plotted are the 
MBH mass, M, , the EMRI mass, m*, the local disk mass, 
jTid, the accreted mass at Eddington rate after lyr (efficiency 
e = 0.1), AM^,AmJ (Sec. IIV A|) . and the accreted mass at 
BHL rate onto the CO for the a and /3-disks per year, Am^^^ 
(Sec. lIVBll . A gap opens outside the radius marked by green 
dots (Sec. lnTB|l . 



Am* are close to the LISA detection uncertainty of the 
EMRI mass, (5to* [Eq. p7|) ] suggesting that the disk grav- 
ity and accretion may lead to detectable effects for LISA 
observations. 

Note that the above mentioned simple formulas de- 
scribing accretion disks are non-relativistic and neglect 
modifications related to the inner boundary condition of 
the accretion disk. If assuming zero torque at the in- 
ner boundary of the disk, ro, this introduces additional 
factors of 1 — (rp/r)^/^ for isothermal disks [3l[, making 
the surface density profile no longer a simple power of 
r. General relativistic corrections introduce additional 
similar factors near the inner most stabl e circular orbit 
(ISCO), light-ring, and horizon [lOl|,[ioi|. Among these, 
the ISCO radius is the outermost one, at 1 < f* < 9 in 
the equatorial plane for spinning BHs. If the shear stress 
generated by the CO heats the disk, then this may fur- 
ther affect the scale height and the density profile. Non- 
axisymmetric inflow across the CO orbit, and the inward 
migration of the CO leads to a more comp licated time- 
dependent density profile [ItI. liol. Hosmoq . We neglect 
these additional factors for simple order-of-magnitude es- 
timates and extrapolate the disk down to r = 3 in many 
of our Figures (which is close to the ISCO for a spinning 
BH with a/A/ - 0.9). 



B. Gap opening 

Up to this point, we neglected the effects of the CO. 
If the CO is sufficiently massive so that its gravitational 
torque moves gas away faster than viscosity can replenish 



it, then an annular gap opens in the disk around the CO. 
The gap width can be obtained from the balance between 
these two competing effects and it is given bj0 



1/3 



(39) 



where 



v^jM, is the dimensionless orbital radius 



of the CO, q = rrii^/Mt is the mass-ratio, and fg is a 
geometrical factor for whi ch STr/ fn « 40 — 50 according 
to numerical simulations |l07l - ll09l |. Typically, the gap 
width is much larger than the horizon radius of the CO. 

Gap opening requires that the equilibrium width f'^ 
be larger than (i) the torque cutoff scale around the CO, 
and (ii) the scale on which the tidal field of the CO domi- 
nates over the SMBH 

[l03,[l0l[Il3[nil- The GO'S tidal 
torque is shifted out of resonance by the mid-plane ra- 
dial pressure gradient and saturates interior to the torque 
cutoff scale, J'cutoff- This is roughl y equ al to the disk scale 
height [see Eq. r^^^^g ~ H [111- The tidal field of 

the CO dominates inside the Hill radius or Roche lobe. 



'Hi 



f'^ - (g/3)i/3f, = 0.32ml/^Af.-^/3r.io. 



(40) 



Gas entering within gets either accreted by the CO or 
it may flow around the CO toward the SMBH. 
Thus, gap opening rcquirej^ 



H < t'a and r'„ < 



Combining Eqs 
the mass-ratio satisfies 



(41) 

we get that a gap opens if 




, (42) 



where in the last line we have utilized Eq. (|30p for v, 
and the terms in the brackets correspond respectively to 



H < 



'A! 



and 



< 



' A- 



Eq. 



p4)) is general for both 
a and ^-disks. For a-disks b — 0, but for /3-disks b = 
1 and the RHS depends implicitly on q and r through 
/3. Substituting H from Eq. (|33]), and (3 from Eq. ((5^ 
and solving for q gives the mass-ratios that lead to gap 



Here, and throughout the paper, primed distances correspond to 
radial distances measured from the CO, i.e. r' = |r — r^-l. 
In the planetary context, the gap opening condition is sometimes 
written as < < where H < guarantees that pres- 
sure effects are less important than the gravity of the CO and 
nonlinearities become significant. Howev er, th e validity and in- 
terpretation of this condition is disputed IllSH , and we shall not 
require it here over the two criteria in Eq. II41I I. 



13 



opening 



I -5/2 . 2 ) 

ggap.c > max<^ OmSa/^^^, 0.092 ai^^ \ , (43) 

9gap,/3 > maxjs.e X 10^* al^^miy^°M^^^^^°f^^^^^° , 
3.9xlO-al/^™^/V/-^-r}, (44) 

for a and /3-disks, respectively, and wc assumed fg = 
0.23. 

Equations (PSI - HI]) can be used to find the CO orbital 
radius at which a gap opens. For a-disks both terms 
decrease quickly with radius, thus the gap may even- 
tually close if the CO orbital radius f^, is sufficiently 
small. However, for /3-disks, the second term depends 
very mildly on radius. Hence, this term for /?-disks is 
best viewed as a radius independent necessary condition 
for the CO mass for gap formation. The CO radius where 
a gap opens, satisfies 



> max ' 



1/2 -1/2 



'''*gap,0 



> 2ial/''mlt^''Mll^''m:,''/'' and 



m.gap,/3 > 3.9 Mq at^'m^/'^M.^f 



(45) 

(46) 
(47) 



Therefore, EMRIs in radiation-pressure dominated a- 
disks in the LISA frequency band (f^ < 50) typically do 
not open gaps unless a < 10^'^. However, a gap typically 
opens around the CO in a radiation-pressure dominated 
/3-disk, provided the CO mass exceeds a value given by 
Eq. (|47p and it is captured by the accretion disk at a large 
radius f^, > 100. As the CO travels inwards, however, it 
eventually crosses the radius [given by Eq. (HH)] where 
the gap starts to refill for /3-disks too. 

If the CO gets first captured in the accretion disk 
around the SMBH at some large radius f^, ^ 100, a 
gap is expected to be cleared quickly, and the gas in- 
terior to the orbit slowly drains down the SMBH on the 
viscous timescale. Depending on disk pa rame ters and 
less understood non-axisymmetric inflows |l03l |. the in- 
ner disk may be completely or partially cleared by the 
time the CO reaches the LISA frequency band at sepa- 
ration fi, < 50. If there is still residual gas interior to 
ri,, the EMRI may eventually catch up with the inner 
disk, shepherding it into the SMBH and causing an EM 
brightening of the AGN [i^. We estimate the radius at 
which this first happens in Eq. (j48p below. Eventually, 
close to the merger, the gap would refill interior to the 
radii given by Eqs. (|45I - H5)) . reigniting the AGN activity. 

The fact that in AGN disks, gaps open around EMRIs 
for large but then eventually close for smaller sepa- 
rations may seem surprising, because it has the oppo- 
site behavior in protoplanctary disks. The reason for the 
difference is the large radiation pressure in AGN disks, 
which makes to be a constant, so that H/r^, decreases 



outwards in Eq. ([55]). In contrast, iJ/r^ is nearly con- 
stant, slowly increasing outwards f or ga s pressure dom- 
inated or self-gravitating disks [H, Il00l | . The other un- 
usual feature in Eq. (|5^ is the factor, where we recall 
P = Pgas/iPgas +Prad)- This factor approaches 1 and be- 
comes unimportant in the gas pressure dominated regime 
for /3-disks where b = 1, while b = makes it identically 
1 for a-disks. However, this factor makes a big difference 
in the radiation-pressure dominated regime for /3-disks, 
where /3 < 1, see Eq. ([211)0 

We note that the above conditions for gap opening, 
based on Eq. (|4T|) arc probably necessary but perhaps in- 
sufficient in realistic disks. While these conditions have 
been well t ested for protoplanctary disks using simula- 
tions |l08l . Il09l |. we are not aware of any studies dis- 
cussing their applicability in 3D radiation-pressure dom- 
inated, turbulent MHD flows for the typical EMRI and 
accretion disk parameters in AGNs. MHD simulations 
of turbulent protostellar disks show that in some cases 
an annular gap may form with an "anti-gap" interior to 
that region where the ga s den sity is increased compared 
to the unperturbed case [ll4| . 



1. Gap decoupling 

As explained above, a gap is expected to form for 13- 
disks in the LISA band for a wide range of parameters. 
The outer edge of the gap at Af can initially follow the 
secondary as long as the GW inspiral rate is smaller than 
the viscous gas inflow rate. As the binary separation 
shrinks, the GW inspiral velocity eventually overtakes 
the viscous inflow rate. Vi,r{fi,) > Wgas. r(Af+) and the gas 
outside the gap cannot keep up with the CO: the outer 
disk decouples Coincidentally, the CO can catch 

up with the disk interior to the gap, if present, causing 
an EM flare [i^. The evolution of the gap and binary 
decouple at < fd, which using Eqs. ([5]) and ([TT]) . is 
given by 

{1.4 X 10"^ a^'^'m~^M~^mliX^ for a-disks, 
15 m!/"a-^/"m.t/"M.-'/"m!/^^A2/i3 
for /3-disks , 

_ (48) 
where in the following we adopt A = 1.7 jlll| . 

The criterion given by Eq. ([48)) is not satisfied any- 
where where a gap has been opened for our nominal set of 
parameters in either a or /3-disks, see Eqs. ([^51 - 115)) . How- 
ever, the gap can decouple in /3-disks for rrii, > 15 Mq or 
a < 0.05 [see Fig. H below] . 



" The 13'' factors were incorrectly missing from the gap opening 
criteria in Ref. [99ll for /3-disks. 
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2.. Density enhancement outside the gap 

As inflowing gas is repelled by the secondary at r^,, 
gas accumulates outside the gap, and the surface den- 
sity is modified relative to its original value without the 
perturber S(f) [Eq. as 



-gap(7-) = S(r) X <( [m^/mdir^] 




if r > Tn 

if fn>f> Tg 

if f < fjr 



(49) 



where fg — Af^ is the o uter boundary of the gajOj a-nd 
k = 3/8 for /3-disk^ [ll5l |. while Vn is the radius at 
which the density enhancement disappears. In practice 
Tn is time dependent; it moves outward with velocity of 
order |wgas,r| |ll5| . By the time the CO enters the LISA 
band Tn > 100. 



In Sec. IIVI through I VIII wc discuss various effects the 
disk has on EMRIs, and show that the opening of a gap 
has a serious impact. Depending on whether a gap is 
opened or not, the CO is subject to Type-II or Type- 
I migration, respectively. Moreover, if a gap is opened, 
then the mass density near the secondary is significantly 
reduced, accretion and hydrodynamic drag effects arc 
quenched. In that case, only the gravitational effects 
play a role (i.e. axisymmetric disk gravity and Type-II 
migration). 



C. Twists, warps, and disk alignment 

In general, MBHs are expected to have a non-negligible 
spin and dominate the angular momentum of the EMRI 
and the disk within r < 10'^. In the absence of an ac- 
cretion disk, radiation-reaction tends to circularize the 
EMRI. Due to GW emission, the orbital angular mo- 
mentum evolves slowly toward anti-alignment with the 
MBH spin, although the total perturbation of the orbital 
inclination is very small [sol. Ill6| . If the CO is on a 
misaligned orbit, collisions with the disk will cause it to 
align or countcralign with the disk [13, 113, . In the ab- 
sence of a CO, if the disk is initially misaligned with the 
MBH spin axis, Lense-Thirring precession will cause the 
disk to warp and twist, and viscous dissipation and ra- 
diative cooling leads to an aligned or anti-aligned config- 
urat i on w ith the MBHs spin (Bardeen-Petterson effect, 
liMIl)- Similarly, a CO can cause warps and twists 



in the disk if on an inclined orbit leading to alignment (or 



Wc assume that the inflow across the gap interior to the CO is 
signiflcantly reduced by the CO. This is a conservative estimate 
since, non-axisymmetric inflow is expected to o ccur across the 
inner edge of the annulus at a reduced rate llOSlI . 
In the radiation-pressure dominated regime, k has not been de- 
termined for a-disks. We conservatively adopt A: = in this 



anti-alignment) of the disk plane with the EMRI jl04| . 
Finally, the accretion disk may be warped by the star 
clust er surrounding the MBH in t he ou tskirts of the disk 
[l20l | , b y an intermediate mass BH [l2l( , or by its own self 
gravity [l22| . Henceforth, we neglect such complications, 
and consider the MBH, EMRI and disk angular momenta 
to be aligned or anti-aligned in the relative radial range 
of LISA observations for simplicity. 



IV. MASS INCREASE VIA GAS ACCRETION 

Next, we concentrate on the mass increase effect due 
to accretion. Our goal is to make order of magnitude es- 
timates and compare them to the detectability measures 
described in Sec. IIIBl 

A. Primary Mass Increase 

The accretion disk feeds matter to the SMBH, but such 
process is bounded by the Eddington limit. This limit 
corresponds to the balance between radiation pre ssure 
and the gravitational force in spherical symmetrjo. The 
corresponding luminosity due to accretion is 



47rGc 



M,, 



(50) 



where M, is the SMBH's mass, which is here the accret- 
ing object, and k is the opacity (see discussion below 
Eq. dMl))- 

If a fraction e of the rest mass energy is converted into 
radiation, then the corresponding accretion rate is 



M, ei\/.c2 



= 2.536 X 10-^efVr"\ 



(51) 



where eo.i = e/0.1 is the normalized efficiency. No- 
tice that the right-hand side of this equation is mass- 
independent. If the SMBH is accreting at a constant 
rate m, = M./M^"^"*, then 



7Tl,l 

ei c' 



(52) 

where A/,_o is the initial SMBH mass. In the analysis of 
Sec, mil wc dropped most factors of ej~^, but as is clear 
from here, every factor of m, is accompanied by a factor 
of ej-^ 

Is such a change in the mass observable via EMRI 
GWs? Equation (|5ip tells us that, during a 1 yr obser- 
vation, the SMBH's mass changes by AM,/M, = 2.2 x 
10~^m,i/ei, which is clearly below the LISA mass mea- 
surement accuracy of Eq. (1171) . Another way to see this 



This limit may sometimes be violated as shown in Sec. IIV Bl 
below. 
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is to compute the phase shift (Eq. for Eddington- 



Umited accretion, which yields 



-Edd 



^^^obs/'/Ec 



10^ m,iM^ ' Tyr , much smaller than the phase 
measurement accuracy. We conclude then that the 
change in the SMBH's mass via Eddington-limited accre- 
tion has a negligible effect on the GW signal, irrespec- 
tive of the type of disk modeled. The accretion has to 
be very supcr-Eddington or the radiation efficiency very 
small (?7i..i/ei — > 10'^) for it to have any impact on the 
GW signal for typical EMRIs at 1 Gpc. 



Eq. ((55)) below). Figure [3] shows the corresponding mass 
accretion rate per year for a and /3-disks, including the 
quenching processes discussed next. 



C. Quenching of BHL Accretion 

Accretion can be "quenched" or suppressed by several 
different astrophysical processes. In this subsection, we 
summarize all such quenching effects that severely modify 
the accretion rates quoted in Eq. ([Si 



B. Secondary Mass Increase 

The CO itself increases in mass too, as it feeds from 
the ambient gas in the accretion disk. We consider the 
case when the CO orbital inclination is aligned with the 
disk. The thickness of the accretion disk [Eq. ([M)) ] is 
much larger than the horizon diameter by a factor of 
iJ/m^ ~ 1.5m,iM,/m* ~ 10^, and thus, it completely 
surrounds the CO. 

In such circumstances, accretion can be ana- 
lyzed within the framework of Bondi-Hoyle-Lyttleton 
(BHL) (29I43H . The characteristic radius of accretion 
can be calculated as the radius at which the thermal en- 
ergy of particles is less than the gravitational potential 
energ-\Pl 



10 



(53) 

The second and third equality correspond to corotating 
quasi-circular COs, neglecting Wroi and using Eqs. (|38p 
and (|34l). 

Assuming isotropic accretion and an adiabatic equa- 
tion of state, the Euler and continuity equations can be 
integrated to give 



\3/2 



1.5 X IQ-"^ yi:-'^ a^^m~^ M~^m^if% 
for a-disks, corotating, circular CO, 



(54) 



3.5 X 10^4 yr- 



-4/5 . -17/5 



M. 



-4/5 



39/10 



for /?-disks, corotating, circular CO. 

where p is the ambient density, Cg is the sound speed, 
and Wrei is the relative velocity of the gas with respect 
to the medium, Eqs. ([8]) and ([36H38| . The numerical 
values shown in the second line are only representative, 
they assume Urei = I'Wgas ^ 'i'*! ^ c-s- This is approx- 
imately satisfied near the SMBH, but in the numerical 
calculations we substitute the estimated value of t;,oi (see 



We use r' to distinguish orbital distances measured from the CO. 
We denote the CO orbital radius by r unless it leads to confusion, 
otherwise . Over-bar denotes units of M, . 



1. Quenching by wind and tidal effects 

When Vrei is not neglected, the estimates of Eq. ([5^ are 
reduced. The relative velocity between the CO and the 
gas contains contributions from the relative bulk motion 
of the gas, differential rotation of the disk, and turbu- 
lence. 

The effect of differential rotation can be estimated as 
follows. Since the gas velocity is different at the edge of 
the accretion range relative to the bulk velocity at r at 
orbital radii r ± Sr' , then 



5vdr 



-3/2 



(55) 

where in the second equality we have set {6r'^,6r''^) ~ 
(r^ , 0) , and in the third used = y^AI,/r and F^^ = 
— 1/r for the Christophel symbol in flat space. Substi- 
tuting in for the Hill radius [Eq. (|40|)]. this becomes 



0.015 



™*i 

M.5 



1/3 



-1/2 
"10 ■ 



(56) 



We estimate the radial wind using Svr ~ 
where Wgas.r and Vi,r can be found in Eqs. 
The relative velocity induced by an azimuthal wind Sv^ 



■^gas.r L'^r | ; 

Bl) and (1371. 



= — Cs 



0.0053 TO^iJ'io^^^ for a-disks 
for /3-disks 



0.013 m^ifi 



' 10 
-5; 
•1' 10 



(57) 

which we derive in Sec. IVl below. 

With all of this in mind, the relative velocity is then 



(58) 



Figure [2] shows that 6vdi dominates the relative gas ve- 
locity relative to the bulk azimuthal and radial wind 
velocities. Indeed, Svdr > Cg when the CO's mass is 
larger than > 300 M^^mlif^Q or if f > fdr = 



31m,iM, 



1/3 -1/3 



In that case, one might expect large 



deviations in Eq. ((54| . As we show in Sec. IIV C 31 how- 
ever, this quenching mechanism is typically superseded 
by limited gas supply. Nevertheless, we include Vrci in 
our calculations below. 
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The accretion is very anisotropic due to differential 
rotation and turbulence in the accretion disk. The accre- 
tion rate in such an advection dominated acc retion flow 
may be significantly less than the BHL rate |l23l . Il24l | . 
We consider our simple estimates to be accurate only to 
the order of magnitude. 



Is Eq. (|63|) an important constraint for EMRI sys- 
tems observable by LISA? If the mass-ratio is extreme 



77 



< 10- 



then 



(Fig. El), so Eq. dMl) im- 



poses a constraint if and only if < M, is violated 
for the radial separations covered by the CO during the 
observation. This is the case outside r > fq, where 



5.1 



2. Quenching by thin disk geometry 

Spherical BHL accretion is valid in the region where 
the Bondi radius is less than (i) the scale height of the 
disk H [Eq. ((Ml)], and (ii) the Hill's radius or Roche 
lobe where tidal effects from the SMBH are negligible 
[Eq. (PU)) ]. These constraints are satisfied interior to 



niiin = 19m,iA/,Y^m^i^^^ , 

or ■ -1/3 
''tidal = 26 m, lA/,^ 



(59) 
(60) 



Thus, the thin disk requirement is always more restrictive 
for radiation-pressure dominated disks. 

Beyond orbital radius fthin [Eq. ((59)) ]. the accretion 
cross section is reduced from Airr'^ to Airr'^H, because of 
cylindrical symmetry. Consequently, the accretion rate 
is modified as 



H 



1, 



3. Quenching by limited gas supply 



(61) 



1/3 -1/3 



8/39^^44/39^6/13 



for a-disks 

-20/39 



for /3- disks . 



(64) 

Figure [2] shows that the accretion rate of the CO sat- 
urates at M, at large f. BHL accretion is supply lim- 
ited for at least part of the observation for the particular 
EMRI systems we consider in Sec. lIXAl below. 



4- Quenching by radiation pressure 

Since the BHL accretion rate is typically super- 
Eddington for a mass nii,, does radiation pressure quench 
such large rates? One has to be careful about this point 
since the derivation of the BHL accretion rate restricts to 
adiabatic flows, neglecting the effects of radiation pres- 
sure, heat transport, and cooling. Super-Eddington mass 
accretion onto the CO is possible if the radiation is trans- 
ported in ward wit h the inflow faster than how it can dif- 
fuse out [12M27I. We compare the diffusion time with 
the infall time below. 

The infall time of the fluid element from the Bondi 
radius to a distance r' from the CO is approximately 



The BHL accretion rate might also be limited by the 
amount of gas supply near the CO. First, note that the 
radial inspiral velocity is typically much slower than 
the radial inflow velocity of the gas Var and VjSr for a and 
/3-disks (Fig.[2|). This implies that there is a constant gas 
flux across the CO orbit from the outer regions. However, 
if the Bondi rate in Eq. ((54)) was greater than the radial 
gas flux toward M, , then the accretion onto would be 
limited by the rate at which gas flows in from the outer 
regions. The mass flux across the CO's orbit is 



27rrS|u„ 



1 



(62) 



where Vi,r and Wgas.r are given by Eqs. ([5]) and ([57]). Thus, 
the accretion rate onto the CO becomes 



Mfl 



UX,*5 ^+ 



(63) 



At large separations, \vi,r \ ^ I'ygas.rl, the CO accretes in- 
flowing gas from the outside and the accretion rate 
becomes independent of radius (as long as a gap is not 
opened, see Sec. lIIIB]) . At very small separations, given 
by Eq. (|48)) . \v^r\ > |wgas,r|, the CO sweeps up the disk 
interior to its orbit, and the accretion rate becomes sen- 
sitive to the assumptions on the available gas supply in- 
terior to the orbit. 



1 - f 



(65) 



where f' = f' I'I''-b is the dimensionless radial distance from 
the CO, and the denominator is the RMS gas velocity. 
The diffusion time from a distance r' from the CO to 



tdm{r') 



(66) 



where we recall that k is the gas opacity (see discussion 
after Eq. ()29|) ). To avoid confusion, wc label the local gas 
density in the close vicinity of the CO as p*{r'), which at 
the Bondi radius is equal to the mean disk density near 
the location of the CO, p*{r'g) = p{ri,) in Eq. ([M)) . For 
BHL accretion, the density increases toward the CO as 
p4r') ^ p^ir'^y-^'"^ [2i-il|. Thus, Eq. ([66]) simplifles 
to 



tdiff = np{r.y^ (1 - f 1/2) 



p;'i/2 



(67) 



where in the second line we have used Eq. ([54| to relate 
tdiff to the BHL accretion rate. Comparing Eqs. ([55)) 
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and (p7)) . we find that the diffusion time is larger than 
the infall time, precisely if > m"'' where 



K m 



B 



?'l/2 



5.6 X lO^^yr^^TO^f^f^o (1 



2-K l+f'l/2 



(68) 



2^:3 



^'1/2 



The second line corresponds to Wrci ^ Cg, but in the 
numerical calculations we substitute Wroi from Eq. (j58p . 

Interestingly, radiation pressure does not have any im- 
pact if the BHL accretion rate exceeds the limit, m"'*. 
For the nominal parameter values for a /3-disk, Eqs. ((54)) 
and ([55)1 show that to^ > m"'* is satisfied for all < ^' < 
1. Therefore, we are reassured that radiation is trapped 
and advectcd inward in this case. However, a-disks are 
much less dense, and this condition is violated interior to 

More generally, BHL accretion may be quenched by the 
various other effects discussed above, modifying r'^ and 
decreasing the gas density, and thereby the diffusion and 
infall times. In this case, radiation pressure may further 
suppress the accretion rate onto the CO if m^, < m"''. 
This criterion can be fulfilled by both a and /3-disks. The 
accretion rate then becomes 



if m'i > mf 
i^'^'^ otherwise . 



(69) 



where m" is given by Eq. ([61]) and we model the 
radiation-pressure quenched accretion as Eddington lim- 
ited, replacing A/, with in Eq. ([5T|). We here choose 
f' = 1, as this gives the most conservative (smallest) es- 
timate for m"'. 



5. Quenching by gap formation 

If the tidal torques of the CO are sufficiently strong to 
dominate over the viscous inflow, an annular gap forms 
around the CO, where the gas density is signiflcantly re- 
duced (sec Sec, nil B|) . Gap formation requires and r^, 
to be sufficiently large, [Eqs. (|45lH71) ]. These conditions 
can be satisfied for /3-disks during the final year of the 
inspiral, but not for typical a-disks. If a gap forms, the 
accretion onto the CO ceases. 

For large CO masses m* > 15 M0 or a < 0.05, the 
inspiral rate becomes faster than the viscous inflow rate 
of gas outside the annular gap if Tgap < < rd (see 
Eg. HS)) . In this case, the CO may "catch up" with the 
gas interior to the orbit [i^ . The inner disk may be filled 
by non-axisymmetric or three dimensional overflow |l03| . 
In fact, in turbulent MHD disks, the region interior to 
the annular gap may have an over-de nsity ("anti-gap") 
relative to the case without an EMRI |ll4l | . In this case, 
rhi, may be restarted interior to r^, and may exceed the 
BHL rate of the original unperturbed surface density of 
the disk [Eq ([M]) ]. However, it is also possible that the 
inner disk drains away before is reached, implying no 
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FIG. 4. Critical CO mass as a function of CO orbital radius 
for various mechanisms to quench BHL accretion onto the 
CO. The accretion rate is reduced for larger or larger 
r*. Top and bottom panels correspond to M, = 10^ and 
W^Mq, respectively. Gap-decoupling occurs interior to the 
green curves. 



accretion. We conservatively assume no accretion onto 
the CO if a gap is present. 



m'l' if r < Tg. 
otherwise . 



(70) 



6. Summary of quenching processes 

The mass increase of the CO is very sensitive to the 
complicated details of accretion disk astrophysics. Most 
of these processes act to decrease the accretion rate from 
. We summarize the EMRI parameters where various 
quenching mechanisms are in play in Figure SI This fig- 
ure depicts the minimum CO mass to^, and orbital radii 
where particular processes become significant to quench 
the BHL accretion rate onto the CO for ai = m,i = 1 
for different M. = 10^ Mq (top panel) and 10^ Mq (bot- 
tom panel). For these parameters, accretion is first 
completely quenched by gap formation for /3-disks, but 
gaps do not form for a-disks for EMRIs in the LISA 
range. Then the gap refills, and accretion is limited by 
the amount of inflowing gas, radiation pressure, differen- 
tial rotation, and the thin disk geometry. Closer to the 
SMBH, these processes become less and less significant 
and the accretion rate increases to the BHL rate m^. 
This decreases inward with the decrease of gas density 
and the increase in the sound speed. The corresponding 
CO mass increase is shown in Figure [5] above. 



D. Implications of BHL Accretion 

Let us now describe the implications of BHL accretion 
on EMRI formation and the GW phase. The former is 
relevant to understand whether EMRIs can remain ex- 
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treme mass-ratio systems as they inspiral in the accre- 
tion disk toward the LISA band, or if they grow in mass 
to form an intermediate mass BH. Then, we estimate the 
corresponding perturbation to the GW phase and discuss 
its detectabihty with LISA. 

Equations ([Ml), (ED), dSSl), and ^ show that 

the CO mass growth rate is sensitive to the location of 
the CO, i.e. is not a constant. The mass at radius f 
can be estimated as 



= m^^o + / mi,dt = m^^o + / , (71) 

J to 



where m*_o is the initial CO mass, and is the accretion 
rate and f is the inspiral rate ([5]). 



1. EM RI formation scenarios 

Does the CO mass grow significantly in the disk prior 
to the LISA observation? If the CO migrated through the 
disk from very large radii, ^ lOM, could be expected 
by the time the CO reaches detectable separations for 
LISA observations 0, [l^l- At large orbital radii the 
accretion rate is supply limited. Assuming that a gap 
does not form and the CO consumes the inflowing gas 
completely and that \ri,\ \vgs,s.r\ 



m M ^ 

Am, < M. iniergcr < 1.2 Mq " f^Qp , (72) 



the starvation mas^J |l29| . which from Eq. 



Am^Q < max 



{ 1800 Mo aj/'m^f M.5f-'o^', 



92OOM0aim2^Af.5fro'} , (73) 

A / fooivj 2/5 . 17/10 ,,r9/10— 29/20 

Am*^ < max 36 Q!;^' to,^' M.g r^Q ' 



(74) 



Note that both Eq. ([72]) and ((721) must be satisfied for 
a-disks, and Eq. ([71)) for /3-disks. In both cases, the CO 
mass remains small ^ 100 for M, = 10^ M©. 
It grows to at most ^ 76 until reaching f = 25 for 
M, ^ 10^ Mq in a /3-disk, but can grow beyond 100 Mq 
for M, > 10^ M0 in an a-disk, if captured in the disk 
outside f > 50. If so, an initial EMRI would morph into 
an intermediate mass-ratio inspiral before entering the 
LISA band. Growth beyond a mass given by Eq. ([73]) is 
halted by gap formation [l7{ . 

We conclude that EMRIs can remain extreme in mass- 
ratio on their journey to the LISA band (f < 50) for 
a non-negligible set of disk parameters. Conversely, the 
mass measurement with LISA could have interesting im- 
plications on the structure of the accretion disk. Suppose 
LISA measures to, to be large, consistent with either 
Eq. (|73p or ([7^ . This information alone would suggest 
that the CO has grown by accretion in an a or /3-disk, 
and suggest the possible presence of a disk, even without 
a direct GW phase shift measurement. In the opposite 
case, if LISA measures to, to be larger than what is ex- 
pected from the growth arguments given by Eq. ([75H7^ . 
then this could point to the common existence of inter- 
mediate mass BHs, which is debated at the time of the 
writing of this manuscript. 



where f2oo = f^o/200 and we have substituted the inspi- 
ral time to coalescence ^merger given by Eq. ([T3)) F^ The 
merger time is typically less than the GW driven inspiral 
time due to angular momentum exchange with the disk 
discussed in Sec. IVIII 17, 35j, so that Eq. (|72p is only an 
upper limit on Ato,. Note that in this case, the mass 
increase is larger than the instantaneous amount of disk 
mass within a few Hill's radii [Eq. ([HI)) ], called the iso- 
lation mass, due to the viscous inflow of gas from the 
outer regions across the orbit. This situation is most rel- 
evant for a-disks where a gap does not form easily [see 
Eq. ([45])]. 

Once the CO mass has grown sufficiently, the tidal 
torque of the CO eventually opens a radial gap in the 
disk, halting further growth. This leads to a limit called 



In the opposite extreme r* 2> Vga.s,r, the CO growth is Hmited 
by the interior disk mass, which is typically less than lOAf^g Mq 
within f ~ 10^, see Fig.|3] 



2. GW Observations 

BHL accretion changes the GW inspiral rate of 
the EMRI due to the increase in the radiating mass 
quadrupole. This leads to a GW phase shift relative to 
a constant to,. Is BHL accretion measurable for EMRIs 
with LISA observations? 

For a crude first estimate let us consider the corre- 
sponding limits on the mass and timescales of Scc. IIIBTI 
Figure [3| shows the BHL mass accretion rate as a function 
of radius (red lines). For these masses (10 M©, 10^ Mq), 



^'^ If a gap opens and the CO is transported inward by Type-II 
migration together with the flow (see Sec. lVIll bclow). then Am, 
is limited by t he lo cal disk mass in a few Hill's radii, i.e. the 
isolation mass [128| . However, typical CO masses exceed the 
local disk mass within f < 10"^ (see Fig. [Sjl, and the inward 
migration rate is slower than the gas inflow rate. In this case, 
gas can build up near the edge of the gap and cause the object 
to grow to the starvation mass. 
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/S.mi,/mi, ^ 10~^ to 10~^ for a and /3-disks. This is com- 
parable to the mass measurement precision of LISA for 
a source at ~ 1 Gpc, suggesting that the perturbation 
caused by the CO mass growth may be marginally sig- 
nificant. The timescale argument with 7^ ^ m^,/m* is 
not directly applicable as to^, varies significantly during 
a one- year LISA measurement. 

A more accurate analytical estimate is computed in 
Appendix \^ where we integrate the total perturbation 
to the GW phase assuming that = Af^, where 
A and B are constants. The phase shift will be pre- 
sented in Scc. lVml Eqs. (fTT5HTT6)) and Tablc|lll it shows 
that the phase shift accumulates with time initially as 
5(j) Ri aiT^ + a2T*, for short observations relative to the 
inspiral rate, and eventually asymptotes to a^T'^"', where 
(ai, a2, 03, C5) are constant coefficients that depend on 
the EMRI and the accretion disk parameters. 

We find that the quenching by gas supply has a major 
effect on the GW phase shift. Without quenching, the 
effect would be of order S(j)^„ ^ 7 and 3000 rads for a 
and /3-disks respectively, and even larger for larger to^,, 
assuming Tobs = 1 yr and M,5 = m^fi = m, = 1. How- 
ever, at most separations, the BHL rate is significantly 
suppressed by gas supply for /3-disks, but the reduced 
dephasing is still about 13rad per a year, and larger for 
larger Af,. The phase shift is proportional to the follow- 
ing combination of accretion disk parameters: ai^m~i 



and 



-4/5 . -17/5 



for the unquenched BHL rate for a 



and /3-disks, and m,i for the supply limited rate. These 
combinations may be marginally measurable by LISA ob- 
servations, given sufficiently strong signals. 



V. HYDRODYNAMIC DRAG 

Next, we consider the drag induced by a difference be- 
tween the gas and CO velocities, sometimes called wind. 
This relative velocity is a consequence of a pressure gra- 
dient in the disk and results in a force that pushes the 
CO both azimuthally and radially [H,!!!]. As in SecHVl 
the goal is to make order of magnitude estimates on the 
corresponding GW phase shift and compare them to the 
detectability measures of Sec. IIIBl 



A. Azimuthal wind 

If the orbital velocity of the gas is different from the 
CO's orbital velocity, the latter will experience an az- 
imuthal headwind or backwind (relative to its unper- 
turbed azimuthal motion). To estimate the orbital veloc- 
ity of the gas. let us write the radial equation of motion 
assuming that M, 3> rUi,, the orbital velocity is much 
larger than the radial velocity, and that the flow is sub- 
sonic: 



M. 



-Vr$^-Vr$disk- 



, (75) 



where and Vr are the azimuthal and radial velocities of 
the gas, Plot is the total pressure of the gas in a comov- 
ing frame, and p is the gas density. For standard thin 
disks, ptat = pel = pH^fl^, where the scale height H 
is independent of radius in the radiation-pressure dom- 
inated regime [see e.g. Eqs. and ([55)) ]. The orbital 
average gravitational potential of the CO, $disk, acting 
on the fluid element at radius r is 



TT r + ?v 



(76) 



where K{k) = /J'^^(l - sin^ Oy^/^dO is the complete 
elliptic integral of the first kind. The gravitational po- 
tential of the disk, $disk, will be given in Sec. I VII [see 
e.g. Eq. dMl)]. 

Let us parameterize the radial density profile via 
p{r) (X r'^ , where the exponent 7 = 3/2 and 7 = —3/5 for 
a and /?-disks respectively [Eq. ([33)) ]. Note that w^VrWr 
is negligible in Eq. ([75)) since Vr ^ Cg and Vr oc l/(rE), 
see Eq. ([57)) . and Fig. [21 Equation ([75)) then becomes 



V 



V.cf>,-V.$,,,, + (3-7):^ 



r 

which one can solve to obtain 



Ah 

r 



1 + (7 - 3)^ + T^V^($* + $disk) 



(77) 



(78) 



In the following we neglect the effects of the potential due 
to the secondary, an approximation valid if the gas ac- 
cretes onto the CO from outside the Hill's sphere, i.e. the 
Tq < 7'^, see Eqs. ([60)) . as well as the disk gravity. 

We then find that a corotating CO always experiences 
an azimuthal headwind with velocity (i.e. orbital velocity 
of gas with respect to the CO), 



5v^ = ilvac?' - 



M, 3-7._ . .2- 
2\ ^^(l-5™.i r 

2 r \ r 2 



(79) 



where we have used Eqs. (1341) and ([781) . This equation 
agrees with Tanaka et. al. [13(j | or the approximate equa- 
tion of Levin [l3|. This estimate, however, does not hold 
for transsonic flows, as in this case the CO experiences a 
backwind, as found by Chakrabarti [2^ . Since such flow 
requires very high accretion rates where the thin disk 
approximation may not hold, we ignore this possibility 
here. 

An azimuthal headwind leads to additional dissipation 
of the CO's specific angular momentum, ^wind = tP /rri-^^ 
where P is the rate of change of the linear momentum, 
so that 



f-wind 



^ _3_7rnJVf. ^^^^^^ 2-3/2 

771^ 2 777,^, 

(80) 

where we have used Eq. ()79p . Clearly, this is typically 
a small perturbation relative even to the loss of angular 
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momentum through GW emission. For unsaturated BHL 
accretion = m^, substituting Eq. ([54]) yields 



wind 



6.1 X 10-^°a^ ^"i.j^M.sffo for a-disks, 



3.4 X lO-<'a-^'''m;;'° M°i''fll' for /3-disks. 

(81) 

For supply limited BHL accretion = M,, Eq. (p^ . 
yields 



/jsup 
wind 



Bondi 



1.0 X IQ-'^ 
2.5 X IQ-'^ 



X TO, 



.3 M^5^2 



' 10 I 



(82) 



where the top and bottom rows correspond to a and /3- 
disks respectively. The change in the angular momentum 
dissipation rate leads to a modified inspiral rate that we 
discuss in Sec. IV Dl 



B. Radial wind 

In addition to the azimuthal headwind, the CO also 
experiences a wind in the radial direction. The corre- 
sponding force is 



-^wind — ^^gas,r 

M,ml 



Hch 



2.5 X m-^^m-{M-^'mij^r[ 



(83) 
(84) 



where we have assumed \vri,\ ^ |wgas,r|, see Fig. [2l The 
radial equation of motion is thus 



2 m^M, 

TO^si r = — 



P. 



wind • 



(85) 



The first term on the RHS is the gravitational force, 
which satisfies m^.M^/r'^ ~ rrii^^ll^^r, by the definition 
of r^vac- Therefore, unlike the azimuthal wind, the ra- 
dial wind does not dissipate angular momentum, but it 
modifies the orbital velocity fivac relative to Keplerian: 



12 q2 

vac 



2.5 X 10-i6!!^ii^ ^ (gg) 



The impact of this modification on the GW phase will be 
discussed in Sec. IV Dl 



C. Dynamical friction 

Dynamical friction is generated by the gravitational 
interaction of a perturber trav eling at some relative ve- 
locity in an ambient medium jl3l| . The gravity of the 
perturber deflects the particles of the medium and gen- 
erates a density wake trailing the perturber. In turn, the 
gravitational pull of the density wake acts like friction, 
decreasing the speed of the perturber. This process is 
analogous to Landau-damping in plasma physics and is 
also important in galactic dy namics fo r objects moving 
through a population of stars |l32l Il33| . 



The standard treatment of dynamical friction in a 
gaseous medium usually assumes that the medium has 
a spatially uniform initial velocity distribution relative 
to the perturber. If the perturber moves on a linear 
trajectory with a subsonic relative velocity in a gaseous 
medium, then the density wake in front and behind the 
pertu rber approximately cancel, leading to a small drag 
force 13J]. However, dynamical friction is more signifi- 
cant for a supersonic perturber. Dynamical friction has 
also been investigated for perturbers moving on quasi- 
circular orbit s in an initially static (i.e. non-rotating) 
medium jl35| . In this case the density wake has a spiral 
structure and the drag force is enhanced. A fully rela- 
tivistic treatment was presented in Refs. (ssl . Il36l |. 

Equations ((23), and jTH]) and Figure H show that 

the relative wind velocity at corotating orbits is typ- 
ically non-rclativistic and subsonic in standard a and 
/3-disks. The arguments mentioned above then imply 
that the standard dynamical friction effect is greatly sup- 
pressed 

The relative gas velocity at different radii outside 
r^, ± however, is supersonic [see Eq. ([55]) ] and thus, 
dynamical friction with respect to the gas in this re- 
gion may be significant. In this case, the velocity of 
the medium is mostly due to differential rotation of the 
accretion disk and not the wind generated by the pres- 
sure gradient effects of Sec. IV Al Differential rotation 
causes the density waves to wind up and standard dy- 
namical friction formulas are not applicable. This regime 
has been well studied in planetary dynamics, which leads 
to the phenomenon called planetary migration. In this 
paper, we explicitly distinguish between such migration, 
described separately in Sec. IVIII below, and standard dy- 
namical friction, which as argued above has a negligible 
effect on LISA EMRIs. 



D. GW Implications 

Are such changes in the orbital dynamics measurable 
with EMRI GWs? Let us first consider the effect of an 
azimuthal wind, for which the inspiral rate is increased 
with the timescale 71. wind = ^/^wind: 

TZ « -4.0 X 10« yr a^ml.M^m^Mo > (87) 
for a-disks and 



« -7.1 X 10^ yr a^/^TOlf M.^^to^^"-/- 



(88) 



for /3-disks. We have here used the fact that the specific 
angular momentum is £ = r^fl where f2 is given in Eq. 

These timescales can now be used to estimate the GW 
dephasing. The effect of the headwind is to change the 



Regular dynamical friction may be significant for standard a and 
/3-disks for retrograde orbits [Sal . 
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orbital angular velocity of the EMRI system, i.e. it in- 
duces a i5f2wind- This change in frequency can be ap- 
proximated as (5f2wind ~ {£wind/i)^St, SO the phasc shift 
is proportional to Taw = ^wind/^- Such a change in or- 
bital frequency induces a change in orbital phase, which 
we can approximate as (50aw ^ ^^^wind^t ^ (5t/7Iw)<^Gw- 
For a-disks we find Scp^w ~ 0-03 rads, while for /3-disks 
we find 6(j)^^ « 33 rads, using a typical set of EMRI and 
disk parameters and and initial separation of fo = 20. 

A more accurate measure of the GW phase shift can 
be obtained by integrating the perturbation to the inspi- 
ral rate due to the hydrodynamic drag. The results of 
a similar calculation was presented in Paper I [H^l for a 
general model, of which the azimuthal wind is a special 
case. We provide the details of the derivation in Ap- 
pendix |B] below. The asymptotic analytical solution is 
presented in Eqs. (|115p - (|116p . where the coefficient and 
exponents are given in Table [III For the default param- 
eters, we get that the perturbation to the GW phase is 
S(j)2w 0.03 rad for a-disks, S(j)^^ ^ 47rad for /3-disks, 
assuming BHL accretion is not saturated for the final 
year of inspiral. Including the limitation by the gas sup- 
ply does limit the process leading to Scp^w 0.15 rad for 
a-disks and Scj)'^^ ^ 0.35 rads for /3-disks. The result of 
the full calculation including all other saturation effects 
is shown in Figures [T] and [5] for different EMRI masses 
and final inspiral radii for a one year observation. We 
conclude then that the azimuthal wind affects the EMRI 
signal only marginally for typical parameters. 

Let us now consider the effect of a radial wind on the 
GW phase. Since the GW phase is proportional to twice 
the orbital phase for quasi-circular orbits, the total GW 
phase shift induced by the radial wind is roughly 

Af Til 

^<^Gw' ^^"^ = 2(f^ - f^vac)T'obs « -t: ' 3 ^ Tohs (89) 

^ 'vac^ ^ 

= 5.0 X 10"^" rads M,^m^im,^f\Q . (90) 

This radial wind is thus completely negligible for the 
LISA measurement unless f > 1000.Fl 



VI. AXISYMMETRIC GRAVITATIONAL 
EFFECTS 

The axisymmetric component of the disk gravity in- 
duces several effects in the orbital evolution of an EMRI: 
it modifies the angular velocity of the orbit and the in- 
spiral rate and induces additional apsidal precession for 
eccentric orbits. We examine these effects here in turn, 
making order of magnitude estimates of the correspond- 
ing GW phase shifts. 



Note that m.i <^ 1 is unrealistic for radiatively efRcicnt, thin 
disk models, but even so, BHL accretion would be quenched by 
the limited gas supply, as explained in Sec. IIV C 31 



A. Accretion Disk Potential 

The gravitational potential of a thin disk may be much 
stronger than the isotropic component of the enclosed 
mass, Mdisk(''^)/?'- This is to be expected since a thin 
ring exerts a much stronger force than a spherical shell 
of the same mass, which can point both in or out for ex- 
terior and interior test particles, respectively. Here we 
estimate the Newtonian gravitational potential of a sta- 
tionary planar disk. 

The total potential of the disk is a superposition of 
the contributions of infinitesimal concentric rings of mass 
dm = 27rrdrl]. Using dimensionless radius variables, 

<I'disk(f) - -M. ^{^r,)^K('^^\ dfo. 

Jr„i„ r + ro \r + roJ 

(91) 

where the surface density I](r) is to be substituted from 
Eq. pT]) or ([33)) . We here used the fact that a circular 
ring of mass dm and radius tq generates a potential at 
r as given in Eq. (|76p . Here rmiu and Tmax a-re the inner 
and outer radii defining the radial extent of the disk, we 
assume fmin ~ and Tmax ^ r^ad in practice [Eq. (|3T|) ] . 

For inclined orbits, with the CO angular momentum 
vector at inclination t relative to the total angular 
momentum vector of the CO and the disk, the potential 
generated by the disk can be expressed more conveniently 
with the Legendre polynomials Pe{x) 

oo 

*disk(f) = -^2^[F2£(0)]2P2£(cOSt) 

xM, i](fo)fo^^dfo, (92) 

where f< = min(f, fg) and f> = max(f, fo) |l37j . The 
integrals in Eq. (j92p can be simply evaluated analytically 
for the particular form of E(f) given by Eq. ([55]). In 
Appendix IC 1[ we carry out this exercise and show that 
in the limit f,„in < ?" < f^ax, 

^•a-disk(0 ~ -TrSaoV^max^^ , (93) 

$^i-disk(r)«27rcoS/^or2/^ (94) 

where Cq = 1.38, and Eq. ([55)1 was used to define the 
dimensionless density scales S^o and S^o as M,E(f) = 
Sqo^"^^^ and T,po^~^^^, for a and /3-disks, respectively. 
Substituting fmax ~ frad from Eq. (|3T|) . we get 

r -1.4 X 10-12«-20/21^.^-^13/21^^22/21_3^ 

$disk(?") I for a-disks , 
m 1.2 X lO-^a^'^'rhli'pfJ'fl/' 

I. for /3-disks , 

(95) 

where we have used <I>(r) = r ^ for the SMBH potential. 

Equation (|95|) shows that the disk potential decreases 
outwards and inwards for a and /3-disks. Thus interest- 
ingly, the disk exerts an outward force on the CO for 
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a-disks. In Appendix I C 1[ we show that this is due to a 
strong quadrupolar field generated by the outskirts of a 
radiation-pressure dominated a-disk. 



B. Change in the orbital frequency 

The orbital angular frequency ^{f) is modified due to 
the axisymmetric gravitational field of the disk from its 
value without the disk fivac given by Eq. Equating 
the centripetal acceleration to the gradient of the gravi- 
tational potential we find 



M, rf$ 



disk 



disk 



J.2 

Therefore the orbital frequency is 



dr 



1 



dr 



(96) 



(97) 



a negligible amount of apsidal precession in the eccentric 
case too. 

The radial oscillation frequency for a nearly circular 
orbit in a perturbed potential is 



fdisk 



disk 



dr^ 



r dr 



(103) 



(see, e.g. jl33l |). li k ^ Q, the ellipse precesses in its plane 
at a rate Q — k. The dimensionless apsidal precession 
rate relative to the Keplerian frequency from Eqs. ((97|) 
and (fTOS]) is thus 



r" I d^$disk , 2 d$disk 



dr'^ 



dr 



(104) 



We relate fiap/i^vac to the induced phase shift below; this 
quantity is proportional to $disk/'& ^ 1- 



E. GW Implications 



C. Change in the inspiral rate 

The binding energy in the Newtonian approximation 

is 

E = -rt'r'7]M -1 + m,$disk , (98) 

2 r 

where M is the total mass and rj = M,7Jii,/AP is the 
symmetric mass-ratio. For EMRIs, r/M « and A/ « 
M,. Assuming no other source of energy loss, other than 
the GW emission, and quasi-circular orbits, 

^^_^^2^,/Vr!6. (99) 
dt 5 ' ^ ' 

We can use Eq. (|97)) to express the radial and time 
derivatives of the energy explicitly as a function of radius 
alone: 

^ = f^fl + 3.-^^+.-^V (100) 
dr 2r^ \ dr dr'' J 

Therefore the radial velocity is modified as 

where rvac is the GW inspiral rate neglecting the effects 
of the disk [Eq. ([5])]. Notice that the first derivative of 
the disk potential does not enter the radial velocity. 

D. Apsidal precession 

While clearly unimportant for quasi-circular EMRIs, 
we briefly show that the axisymmetric disk gravity causes 



The disk potential changes the GW frequency and the 
inspiral rate at fixed orbital radii. The dimensionless 
change in the frequency 5f2 = (fi — fivac)/f^vac, inspiral 
rate Sf = {f — fvac)/?^vac, and the dimensionless apsidal 
precession rate f2ap/f^vac are all proportional to deriva- 
tives of <I'disk('^)/*I'(^)i up to a factor of order unity [see 
Eqs. (|97ll02|104p ]. Equation ^ shows that this ra- 
tio is typically very small, at the level of $disk/*I' ^ 
(2 X 10~^, 10^^) for (a, /3)-disks with very small a and 
large M, and f, {a,M,,f) ~ (10"^, 10*^ M©, 100), and 
even smaller for more typical values. This leads to a small 
change in the GW phase of Eq. (HH). In Appendix [Cl we 
carry out a detailed calculation that shows that for both 
a and /3-disks, self-gravity effects induce a dephasing of 
approximately 10~^ — 10~''rads for our nominal set of 
parameters in a one year observation, which is clearly 
below LISA'S observational threshold. 

These estimates are modified if the gravitational 
torques from the CO quenches gas inflow onto the SMBH, 
and clears a gap (Sec. IIIIBp . In this case, the gas den- 
sity is greatly reduced interior to the CO, but extra gas 
accumulates outside its orbit, changing the potential. In 
Appendix [C] we show that this effect increases the de- 
phasing by roughly an order of magnitude for /3-disks, 
but such an increase is still well below the threshold of 
detectability. Axisymmetric gravitational effects may be 
more important for eccentric orbits with large semi-major 
axis f - 1000. 

VII. MIGRATION 

Let us now examine the role of non-axisymmetric grav- 
itational effects induced by the disk, which lead to a phe- 
nomenon known as "migration" in planetary dynamics. 
As before, the goal is to make order of magnitude esti- 
mates of the migration effect on the GW phase and com- 
pare them with the detectability measures of Sec. IIIBI 
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A. General Properties 

Consider the non-axisymmetric gravitational effects of 
the disk, leading to angular mome ntum diss ipation, anal- 
ogous to planetary migration (see |ll3l . ll28j for reviews). 
The orbiting CO exerts a nonzero average gravitational 
torque on the gaseous disk, creating a spiral density 
wave. The total angular momentum budget is dissipated 
through viscosity and the outward angular momentum 
transport of the spiral density wave. The gravitational 
torque of the spiral density wave exchanges angular mo- 
mentum resonantly with the CO causing it to migrate 
jllO| . This effect is analogous to dynamical friction (see 
Sec. IV C]) . but accounts for the inhomogeneous velocity 
field of the gaseous medium. 

In planetary dynamics, this phenomenon in different 
regimes is called Type-I, Type-II, and Type-Ill migra- 
tion. The distinction is whether a gap is opened (Type- 
II) or not (Type-I) and whether strong corotation torques 
related to horseshoe orbits are taken into account (Type- 
Ill). In the following we neglect Type-Ill migration be- 
cause it is most relevant only if the disk mass near coro- 
tating orbits is comparable or larger than the secondary 
mass. Figure [3] shows that this is not the case for EMRIs 
in the LISA frequency band. 

In a pioneer paper, Goldreich and Tremaine presented 
the first study on Type I migrat ion, using 2 dimensional 
(2D) linear perturbation theory [llO[. Their results were 
later improved by Tanaka et. al. to account for corota- 
tion resonances and 3D effects in isothermal disks |l30l |. 
which are also consistent with 3D hydrodynamic simu- 
lations of laminar disks jl38l | . Fur ther generalizations 
exist for locally adiabatic 2D disks |139| . However, re- 
cent MHD shearing box simulations of turbulent proto- 
planetary disks with a low mass planet show that the 
torques exhibit stochastic fluctuations for Type I migra- 
tion, where even the sign of the torque changes rapidly 
[Ha [mi- Stochastic migration is not expected if the 
satellite is more massive relative to the disk. The mi- 
gration is also very sensi tive to the value of opacity and 
radiation processes (l42j . Recently, Hirata generalized 
the original study of Goldreich and Tremaine for laminar 
2D relativistic disks and found the angular momentum 
tran sport to b e larger by a factor ^ 4 close to the central 
BH |l43l Il44l | . Finally, higher order PN corrections may 
resonantly excite persistent spiral density waves close to 
the SMBH even without an m^,, which could modify the 
torque estimates for EMRIs |l45| . 

To our knowledge, Type-I migration is unex- 
plored for the expected environments for LISA EMRI 
sources: for radiation-pressure-dominated, optically- 
thick, geometrically-thin, relativistic, magnetized and 
turbulent disks, where the mass of the perturbing body 
exceeds the disk's mass. For a simple estimate below, 
we explore the GW phase shift in an EMRI system due 
to Type I migration, using the isothermal non-relativistic 
laminar formulas of Tanaka et. al |13C| . where we include 
the effects of radiation pressure by using the correspond- 



ing expressions for the sound speed and the scale height. 

If the CO is sufficiently massive and/or far out, a cen- 
tral region is cleared out and a circular gap develops (see 
Sec. nil B|) . Once a gap opens, the angular momentum ex- 
change becomes much more regular than Type-I and it is 
determined by the angular moment um t r ansport t hrough 
the more distant part of the disk [HI [l40, \u^. This 



type of angular momentum exchange is referred to as 
Type II migration. If the local mass of the disk is greater 
than the mass of the perturbing object, then the object 
migrates inward on a viscous timescale with the velocity 
of the accreting gas. However, for LISA EMRI sources, 
the local disk mass is smaller than the CO's (see Fig.[3|). 
In this case, the migration slows down, and b ecomes "sec- 
ondary dominated Type-II migration" [l28l |. 

Simple steady-state estimates, based on angular 
momentum balance, were presented by Syer and 
Clarke jll5l |. These estimates account for the increase 
of the gas density relative to an isolated disk, assum- 
ing that the angular momentum exchange is dominated 
by th at ne ar the inner-edge of the disk. Later, Ivanov 
et. al. |104l | relaxed the steady-state assumption and esti- 
mated the quasi-stationary, time-dependent evolution of 
the disk and satellite using a zero stress boundary con- 
dition at the location of the binary. 

Both of these studies focused on thin, one-zone, gas 
pressure dominated, Shakura-Sunyaev disks, where the 
density is an increasing function of stress. This assump- 
tion is satisfied for /3-disks but not for radiation-pressure 
dominated a-disks. In the later case, we are not aware of 
literature that is applicable to Type II migration. Recent 
studies have shown that migration with an annular gap 
is affected by global edge-modes for massive disks, and 
can cause stochastic migration of planets either inward 
or outward (similar to the Type I case with MHD) jl46l | , 
and the migration rate is also sensitive to vortex forming 
instabilities |l47l |. However, these phenomena have not 
yet been explored in AGN accretion disks for EMRIs. 
For an order of magnitude analysis of Type II migration, 
we use both Syer and Clarke and the asymptotic Ivanov 
et. al. equations to estimate the corresponding GW phase 
shifts for EMRIs in /3-disks. We make a conservative es- 
timate for radiation-pressure dominated a-disks, based 
on angular momentum balance between rrii, and the lo- 
cal disk, neglecting the accumulation of gas outside the 
gap. 



B. Type-I migration 

For an isothermal 3D disk. Type I migration removes 
or increases angular momentum at a rate 



±ci- 



uii, r 



(105) 



where ci = (1.4 — O.57) and 7 = 3/2 or —3/5 for a 
and /3-disks, respectively, and the ± signs highlig ht th e 
stochastic nature of migration in a turbulent disk jl30| . 
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The magnitude of ci is different for locally adiabatic 2D 
disks by a factor ~ 2, b ut th e scaling with other param- 
eters remains the same [139| . Relative to the GW rate of 
angular momentum loss, 



^mig,I 



±7.2 X 10 
±5.1 X 10^ 



a 



-4/5 



^M.srfg for a-disks , 



for /3-disks . 



-7/5,. 6/5 59/10 
• 1 ^".5 '10 



(106) 



Notice that the Type I migration dominates over even 
the GW loss of angular momentum at sufficiently large 
radii, beyond f « 150 for a-disks and f « 112 for /?-disks. 



C. Type-II migration 

If the CO is sufficiently massive and/or far out that a 
gap opens (see Sec. IIIIB[) the CO is subject to Type II 
migration. The angular momentum exchange depends on 
the local disk mass near the binary, rrid given by Eq. psp . 
Typically, md < m^, (sec Fig. [3]), and the specific angu- 
lar momentum dissipation rate, in t he quasi-stationary 
approximation of Syer & Clarke jll5| , is 



^mig,II,SC — ( — - 1 ^gas — 



^orb^gas.r 



where €gas is the angular momentum loss in the gas due 
to viscosity, Wgas.r = 2M,r/md(r) is the radial velocity of 
gas given by Eq. ((57)) . Uorb is the orbital velocity Eq. ([5]), 
and fc = 3/8 for electron-scattering opacity and /3-disks. 
If neglecting the banking up of gas near the outer edge of 
the gap, then angular momentum balance implies fc = 1, 
which we ad opt conservatively for a-disks [see eq. (42) 
in Ref. jl04l |]. We note that this assumption is different 
than those used in Refs. [s^ H^. Substituting Eqs. ([6|), 
dnS]), and dSZl), 



Kmig,II,SC 



6.2 X 10-6 

for a-disks 
5.8 X 10 



3 1/2 . 5/8 , ,rl3/8 -11/8-25/8 

a/ TO.i M.J m-i, rir, 



for /3-disks . 

(108) 

Notice that for a-disks, this is independent of a; the dis- 
sipation of the CO's angular momentum is proportional 
to M,. For /3-disks, Eq. (|108p also accounts for the accu- 
mulation of mass near the CO, which leads to additional 
angular momentum dissipation, sensitive to a. 

We also consider the time-dependent solution of Ivanov 
et. al. |l04| for /3-disks that accounts for the accumulation 
of gas and a zero stress boundary condition at the edge 
of the gap (see also Eq. (24) in Ref. |l48l |): 



•^mig.IIJPP 



-C2«orb — (r rdd)'/', (109) 



where rdd is the "radius of disk dominance" , which sat- 
isfies md(rdd) = [see Eq. ([551) ]. C2 = {1 + S[l — 



(r/rdd)^/^]}''", S = 6.1, and fca = 0.26. Note that C2 
is only mildly r dependent. For r ~ rdd, C2 = 1 so 
that both estimates [Eqs. (|107p and (|109p ] imply a mi- 
gration rate tracking the radial velocity of gas, so that 

•^migJIJPP ~ ^mig,II,SC ~ ^gas = " I'orb'^gas.r/S, whilc for 

r ^ rdd, C2 ~ 1.66, asymptotically independent of r. 
Similar to Eq. (|107p . Eq. (|109p is valid for electron scat- 
tering opacity /3-disks, but is not applicable for radiation- 
pressure dominated a-disks. Substituting Eq. ([55]) for 
Tdd, we find for r < rdd, 



' 10 



(110) 



D. Quenching of migration 



In the following we assume that inside and outside the 
minimum gap opening radius Tgap [Eq. (^SHiS)) ] Type-I 
and Typc-II migration operate, where we use Eq. (jl06p 
([TU5HTTU)). respectively. Once the CO has 
Tgap, the gas is no longer expelled effi- 



and Eqs. 
crossed inside fg 
ciently by the tidal field of the CO, and the gas is free to 
flow in on the accretion timescale with a radial velocity 



, (107) V: 



[Eq. ([57])]. If this is faster than the GW inspiral 
Vi,r [Eq. dH)], then the gap can refill. In 



gas,r 

rate of the CO 



Sec, nil BT] we have shown that \v^r\ < ["'^gas,? ! is satis- 



gas,r| 

fied if f > fd, given by Eq. (|48| . Thus, if f = fgap > fd 
is met [see Eqs. (|^51 - H5)) ]. we assume that the gap refills, 
and switch from the disk generated torque from Type-II 
to Type-I. Otherwise, if a gap has formed and cannot fol- 
low the inspiral rate of the CO, disk torques are expected 
to shift out of resonance, and become too distant for ef- 
ficient angular momentum exchange. Then the interac- 
tion is greatly suppressed, and we assume ^mig(?^) = 
if ^gup ^ r < fd- Once the gap has formed and de- 
coupled, it can no longer refill, so we further assume 
^mig = if f < 7"gap < fd. Finally, if the gap refills 
(i.e. fgap > fd) then Typc-I migration can operate effi- 
ciently even if the inspiral rate is faster than the viscous 
inflow rate (f < fd), since the spiral patterns can form 



much faster on a dynamical timescale 



dyn 



n- 



The 



interaction may be significantly different only much later 
when the inspiral time tew = f/f aw becomes faster than 
the cooling time of the disk Il49l . i.e. 



cool 



(111) 



where few is the GW driven inspiral rate, fi is the orbital 
angular velocity and a is the parameter in the viscosity 
prescription. From Eqs. ((2)), (E)), and (|llip we get few < 
tcooi inside 



0.055 



2/3,,r-2/3^ 



2/3 



(112) 



Typically the disk cooling docs not impose a limitation 
outside the ISCO for EMRIs unless a is very small and 
the mass-ratio is not extreme. 
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To summarize we assume, 



J ^mig,I J if r < Tgap , 



'-mig,!! : 



' gap 1 



and 



i-mig 



where for 



i' ■ 

mig 



if r < rc or '"gap <. r -C 
or f < fgap < rd , 
otherwise , 



(113) 



(114) 



^mig,i wc use Eq. (|106p . while for ^migji ei- 
ther the Syer-Clarke model Eq. (|108p or the Ivanov 
et. al. model (|110p . The later is only available for /3- 
disks in the radiation-pressure dominated regime, while 
the Syer-Clarke is applicable both for a or /3-disks. In 
summary, we consider three cases for migration, which 
utilize the same Type-I model, but differ in the Type-II 
regime 



1. M asc: a-disks with Syer-Clarke Eq. ([TU5)) : 

2. M /3sc: /3-disks with Syer-Clarke Eq. ([TU5)) : 

3. M /3ipp: ^-disks with Ivanov et. al. Eq. ([TU5)) . 
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TABLE I. Summary of accretion disk effects on the GW phase 
shift induced by different accretion disk effects relative to 
vacuum waveforms. Rows correspond to different accretion 
disk effects: Primary Eddington-limited accretion (Primary 
acc. ), secondary Bondi-Hoyle-Lyttleton (BHL) accretion, az- 
imuthal winds (W), disk's self-gravity (SG) and migration 
(M). Columns correspond to different EMRI systems assum- 
ing a 1 year observation. The entries x (y) denote x x 10^ 
in radians. The phase shift is negative for all effects except 
for SG a. Observe that /9-disk migration is the dominant ef- 
fect for the first two columns, while /3-disk BHL accretion is 
dominant for the last two columns. 



E. GW Implications 

The change in the angular momentum dissipation rate 
modifies the GW driven inspiral rate, which modifies 
the GW phase evolution. The corresponding phase shift 
(5(/)GW can be calculated in a similar way as for an az- 
imuthal wind (see Sec. lVDI above. and Paper I [H^]). The 
interested reader can find the details of the derivation in 
Appendix [B] 

For the default parameters, we find that Type-I mi- 
gration produces a typical GW dephasing on the order 
of 10"^^ and 10 radians during the final year of observa- 
tion in a and /3-disks, respectively. The large increase in 
dephasing for the /3-disk case is due to the much larger 
surface density. 

Even more interesting is the case where a gap opens; 
Type-II migration is clearly the most dominant pertur- 
bation among all the ones considered in this paper. How- 
ever, as showed in Sec. IIIIBl a gap is expected to close 
for most EMRIs in the LISA band. For a-disks, Type-II 
migration can generate GW phase shifts of order 10—100 
radians in one year, and up to lO'^ — lO'' rad for /3-disks. 
The modifications to the GW spectrum are so large that 
if this effect is in play, that vacuum EMRI templates 
might be ineffective to extract LISA GWs. The blue lines 
in Figure [T] and [5] show the phase shift due to migration 
as a function of final orbital radius including quench- 
ing effects [Eq. (|114p ] for both Type-II migration models 
(Syer & Clarke steady state and the Ivanov et. al. quasi- 
stationary models, Eq. (|108p and (|110p . respectively), as- 
suming ai = m,i — 1. For other system parameters and 
observation times, we provide asymptotic analytical ex- 
pressions in Eqs. (|115H116|) and Table HIl below. 



VIII. COMPARISON OF ACCRETION DISK 
EFFECTS 

We summarize the GW phase shift generated by the 
dominant accretion disk effects in Table ID Different rows 
show different accretion disk effects with the parameters 
of Eq. ((26)) . while columns show different EMRI systems 
with component masses given by all four combinations of 
M, = (10^ 10^) Mq and to* = (10, 100) Mq. The entries 
represent the GW phase shift between the standard vac- 
uum waveform and those including the effects of the ac- 
cretion disk in the Newtonian approximation. The phase 
shift is between one-year long GW waveforms with the 
same final radius in the most sensitive LISA frequency 
bandEB (r/,ro) « (16,25) for (M.,to*) = (10^ 1O)M0; 
(f/,fo) « (16,42) for (Af.,TOO = {10^10^)Mq; 
(f/,fo) « (3,7.6) for (Af.,TOO = (1OMO)M0; (f/,fo) « 
(3,13) for (M.,TO*) = (lO^lO2)Af0. The phase shift 
estimates are derived in Appendix |3 |Bl and [CI using 
Eqs. (jAlOp . ((B8)l . and jCS)) . where the underlying quan- 
tities TO*, £ and $disk are substituted from Sec. llVflVIl 
taking into account all of the quenching mechanisms that 
are in play. In particular, for migration we use Type I 
and Type II in the appropriate radial ranges without and 
with gaps, respectively, where we utilize either the Syer- 
Clarke (SC) or the Ivanov et. al. (IPP) model for Type 



Here we don't marginalize over an arbitrary phase shift between 
the two waveform templates for simplicity. These estimates are 
in good agreement with the more complicated calculations pre- 
sented in Sec. 1X1 below. 
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FIG. 5. The GW phase shift as a function of final radius in units of M, induced by different accretion disk effects relative 
to vacuum waveforms. As in Fig. [1] solid curves correspond to a-disks, while dotted ones to /3-disks. Different color curves 
represent different accretion disk effects (ai = m.i = 1): black is for secondary Bondi-Hoyle-Lyttleton (BHL) accretion, green 
is for azimuthal winds (W) and blue is for migration (M), with dotted and dot-dashed curves for /3-disks in the Syer-Clarke 
and Ivanov et. al. model respectively. The thin, solid magenta line is the total accumulated GW phase in vacuum, while the 
thick, solid (dashed) magenta line is a measure of the sensitivity to which LISA can measure the GW phase for a source at 1 
Gpc (10 Mpc). This figure presents similar conclusions to those found in Fig. [T] but for a different set of EMRI parameters. 



II migration, see Sec. IVII Dp^ 

Table U shows that the migration estimate incorporat- 
ing the SC model in /3-disks is by far dominant, followed 
by migration using IPP in /3-disks, and BHL accretion. 
The effects of the disk's self-gravity, wind effects, mi- 
gration in a-disks and Eddington-limited accretion onto 
the SMBH are ah completely negligible for LISA EMRIs. 
All effects studied lead to a reduction of phase cycles in 
a fixed observation time for fixed final radius, except for 
the effect of disk gravity in a-disks. Other detectable 
EMRIs not shown in this table (different masses and/or 
orbital radii) might be more sensitive to disk effects. 

Figures [T] and [S] show GW phase difference induced by 
various disk effects as a function of final orbital radius. 
Solid curves correspond to a-disks, and dotted curves to 
/3-disks. Curve colors represent different accretion disk 
effects: black curves correspond to BHL accretion, green 
curves to azimuthal wind effects, and blue curves to mi- 
gration (dotted for SC-^, dot-dashed for IPP-^). Left 
and right panels correspond to different mass ratios and 
SMBH masses. Additionally, Figs. [1] and [5] include three 
disk-independent reference curves. The thin, solid ma- 
genta line is the total accumulated GW phase without 
disk effects. The thick magenta curves represent the ap- 
proximate accuracy to which LISA can measure the GW 
phase as a function of radius, using the simple estimates 
given by Eq. (fT9|) for a typical source located at 1 Gpc 
(solid) and 10 Mpc (dashed). We present more detailed 
estimates of measurement accuracy in Sec. 1x1 below. 



Typically migration in /3-disks generates the dominant 
phase shift for all final radii. The large drop in the right 
panel occurs because the gap closes around f ~ 20.4 for 
/3-disks for these masses, changing Type-II into Type-I 
migratiorF^. BHL accretion follows migration in impor- 
tance, and in certain cases (as in the right panel of the 
figure), the former can be the most important effect if the 
gap refills. The sudden quenching features in the BHL 
accretion curves are due to EM radiation pressure, where 
photon diffusion becomes sufficiently short for EM radi- 
ation to escape the flow. Observe that LISA is sensitive 
to EMRIs that are close to the SMBH (e.g. r < 50); the 
sharp rise of the LISA accuracy estimate is because EM- 
RIs at larger radii are not detected with SNR of 10 or 
bigger. All disk effects above the thick magenta line are 
significant and possibly measurable, while wind effects 
and self-gravity effects (not shown in the figure) always 
lead to dephasings of order 1 radian or much less. These 
figures suggest that GW observations may be used to test 
the predictions of different accretion disk models. 

We note that for the masses in the left panel of Fig. [5l 
the CO is located at = 43 and = 51 one and two 
years prior to merger. A gap is open beyond f-i, > fgap = 
96 and 5 for a and /3-disks. The EMRI evolution becomes 
much faster than the viscous inflow and the gap decouples 
at a radius of fd = 39 for a /3-disk. 

To describe the phase shift for different observation 
times or accretion disk parameters than those shown in 
Figures [T] and [5l we provide analytical expressions for the 
asymptotic phase shift for various processes in two limits. 



To understand which quenching mechanisms are in effect for the 

particular cases, see Figure|4]in Sec. IIV C 61 Gap decoupling does not occur since ~ 19.2. 
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TABLE II. Constant coefficients and exponents in Eqs. (|115H116|l for long and short observations (columns), for different effects 
(rows). We use the notation x {y) = x x 10^ in radians. Observe again that tire dominant dephasing is due to Type-II migration 
(M2) both in tire short- and long- observation limits. Also notice that all dephasing scales with positive powers of the radius 
(i.e. C6 > and da > 0), implying accretion disk effects become stronger for EMRIs orbiting at larger separations. 



IX. RELATIVISTIC WAVEFORMS AND 
DETECTION 

Next, we consider accretion disk effects in more realis- 
tic waveform models. EMRI GWs are highly relativistic, 
with velocities close to the speed of light and sometimes 
skimming the SMBH horizon. As such, Newtonian wave- 
form estimates for GW data analysis are inaccurate. Here 
we investigate the relativistic correction to the EMRI 
dynamics using the extended one body framework, and 
make simple estimates on the imprint of accretion disk 
effects on the GW waveform. This analysis, however, 
continues to neglect relativistic corrections to accretion 
disk effects. 



A. Systems Investigated 

In the rest of this paper, we restrict our investigations 
to the following two representative EMRI systems: 

• System I: Masses M, = 10^ Mq, = 10 Mq, 

spin parameter a,/M, = 0.9, observation time T = 
lyr, range of orbital radius f S (16,25), orbital 
velocity v/c G (0.2,0.25), GW frequency /gw € 
(0.005,0.01) Hz, GW phase cpcw 1-3 x 10^ rad. 

• System II: Masses M, = 10^ Mq, to^ = 10 A'/©, 
spin parameter a,/M, = 0.9, observation time 
T = lyr, range of orbital radius f S (3,7), or- 
bital velocity v/c £ (0.37,0.54), GW frequency 
/gw e (0.003, 0.01) Hz, GW phase 0gw 9 x lO^ 
rad. The ISCO is located at nsco ~ 2.32. 

Figures [T] and [S] shows that accretion disk effects are ex- 
pected to be significant for these systems. 



namely 



(115) 



where the (C,c,D,d) parameters are given by the rows 
of Table HIl Here di = c\ and c?2 = C2. The first (sec- 
ond) formula is applicable if the observation time is much 
shorter (longer) than the inspiral time [Eq. ((T^ ]. The 
particular processes are represented by rows in order: 
the full BHL rate (i.e. assuming no quenching throughout 
the observation), the supply limited BHL rate (SB, using 
Eq. inn throughout the observation) , the corresponding 
hydrodynamic drag from the azimuthal wind for the full 
BHL accretion and SB (W and SW, respectively) , Type-I 
migration (i.e. assuming a gap is not present throughout 
the observation), steady-state Type-II migration, quasi- 
stationary Type-II migration (i.e. assuming a gap is open 
throughout the observation), and disk self-gravity (SG). 
Note that LISA observations are sensitive to the combi- 
nations a1^rn'^\, see corresponding columns. 

Most accretion disk effects are several orders of mag- 
nitude larger for /3-disks relative to a-disks. This is be- 
cause /3-disks can be much more massive in the regime 
of interest for EMRIs. In particular, this suggests the 
GW phase shift may be used to test the predictions of 
different accretion disk models. 
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We make the foUowing simphfying assumptions. First, 
we consider only quasi-circular EMRIs on the equatorial 
plane, such that the orbital angular momentum is per- 
pendicular to the SMBH's spin angular momentum. We 
have also investigated Systems with spin anti-aligned or 
zero and found similar results. The accretion disk is also 
assumed to be on the equatorial plane, such that the 
EMRI is completely embedded in the disk. We ignore 
the CO's spin angular momentum, as well as sub- leading 
mass-ratio terms in the radiation-reaction fluxes and in 
the Hamiltonian. 

Such simplifying assumptions make the problem ana- 
lytically tractable within the EOB framework, employed 
here for waveform modeling. As of the writing of this 
paper, the EOB framework for EMRIs has not been 
sufficiently developed for non-equatorial orbits; it has, 
however, been satisfactorily tested for equatorial orbits 
with extreme or comparable mass-ratios, q < 10^'* and 

q > 10-^ 

We expect that many accretion disk effects (migra- 
tion, BHL, and wind effects) will be maximal for the 
equatorial EMRIs studied here. Other effects, however, 
are substantially different for EMRIs inclined with re- 
spect to the accretion disk. Similarly, disk effects may 
excite eccentricity, which may dramatically increase the 
impact of accretion disk effects on the GW observables 
[ia I47I [53I [S^ l . A study of non-equatorial or eccentric 
EMRIs with an accretion disk is beyond the scope of this 
paper. 



The binary's binding energy in a Kerr back ground for 
a quasi-circular orbit is simply given by |l69| 



1 - 2f- 



-3/2 



(118) 



where fj, = Af,m^,/M is the reduced mass, with M = 
M, -I- m* the total mass, f = r/M, and x, = a,/M, is 

the reduced Kerr spin parameter. Notice that it is the 
binding energy that drives the orbital evolution, and not 
the total energy of the system, which would also account 
for the rest-mass energy. 

We employ here the fact orized form of the GW flux, 
considered in jl58l II6II Il62{ , with the assumption of adi- 
abaticity: 
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i=2 m=0 



where R is the distance to the observer, m is the az- 
imuthal quantum number of the multipolar-decomposed 
waveform, and 



him — h'l^''"" Tim e"'*" (pim) 



(120) 



and where Cp is the multipolar waveform parity (i.e. ep = 
if £ -I- TO is even, ep = 1 if £ -|- m is odd). All the 
terms in Eq. i^2Ql (S'/m-- Sim and p t m. ) are fun ctions 
of {r,(f>,n) that can be found in [Hi, |l6l!, [TH]. The 
Newtonian part of the waveform is given by 



B. Basics of the EOB Framework 

We employ the adiabatic EOB framework of [73 
to model EMRI waveform s. The EOB scheme was 
first proposed in (l50l Il5l| to model the coalescence of 
comparable-mass BH binaries. Since then, this scheme 
has been gr eatly enh anced and extended to other type 
of systems [T52l - [l63| . Waveforms constructed in this 
way have been succes sfully com pared to a set of numer- 
•4- 16_6|l a. nd to self-force calcula- 
proposed the combina- 
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Recently, 



73H2E 

tion of EOB and BH perturbation theory techniques to 
model EMRI waveforms for LISA data-analysis purposes. 
This is the scheme we adopt in this paper. 

In the adiabatic, EOB framework for quasi-circular 
EMRIs, the GW phase can be obtained by solving the 
adiabatic equation 



dE 



(117) 



where = is the orbital angular frequency (Eq. ([2]) 
but with relativistic corrections), overhead dots stand for 
time derivatives, E is the binary system's total binding 
energy and J"gw = dEoyj/dt is the GW energy fiux. We 
have here implicitly assumed a balance law: all loss of 
gravitational binding energy is removed only by GW ra- 
diation, dE/dt = —dEcvj/dt. 



Km = 



Ye. 



,{^/2,<p). (121) 



where Yi^miO^'t') arc the standard spherical harmonics, 

n^g^ and Q+£p are n umerical coefficients that depend on 
the mass-ratio |l6l| . The orbital velocity v is related 
to the orbital frequency via v = (Mil) 1/3, which then 
implies the binary orbital separation is 



1 2/3 



^2/3 



(122) 



where recall that the overhead bar stands for normaliza- 
tion with respect to total mass: f2 = MVl « M,Vl. 

At this stage one might be slightly confused, as the 
right-hand side of the evolution equation one wishes to 
solve (Eq. (|117p ) depends on a variety of quantities, in- 
cluding the orbital separation and the orbital phase. The 
assumption of adiabatic quasi-circularity allows us to re- 
place the orbital separation in terms of the orbital fre- 
quency via Eq. (jl22p . By definition, the orbital phase 
is related to its frequency via the differential equation 
(j) = Q.. This equation, together with Eq. (|117p forms a 
closed system of coupled, first-order partial differential 
equations that can be consistently solved. 

The flux in Eq. (jll9p . however, is not sufficiently accu- 
rate to model EMRIs. First, it neglects the loss of energy 
due to the absorption of GWs by the MBH. Second, it 
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is built from a PN expansion, which is in principle valid 
only for slowly-moving sources, which EMRIs are not. 
This flux can be improved by linearly adding BH ab- 
sorption terms and by adding calibration coefficients to 
Eq. (|119p that are fitted to a more accurate, numerical 
flux. This is the procedure proposed in |73l475j , which we 
follow here. We include up to 8 calibration coefficients, 
obtained by fitting to a more accurate Teukolsky evolu- 
tion in the point-particle limit, as given in Eqs. (26)-(29) 
of 75|, as well as BH absorption terms as given in Ap- 
pendix B of [Zl. 

Initial data for the evolution of the system of differ- 
ential equations is obtained through a mock evolution, 
started at r = lOOM, and ended at few = O.Of Hz (see 
e.g. [73l - [75j ). The mo ck ev olution is initialized with the 
post-circular data of |l5l| . Once the evolution termi- 
nates, one can read initial data one- year prior to that 
point directly from the numerical evolution of the mock 
simulation. We obtain initial data of such form sepa- 
rately both in the case of vacuum EMRIs and for EM- 
RIs in an accretion disk, as the evolutions are generically 
different. Once the orbital phase is obtained by solving 
Eq. (|II7p with this initial data, the waveforms are readily 
obtained through Eq. ([1^ . 

Before proceeding, let us comment on the differential 
system one must solve numerically. As already men- 
tioned, since the source of Eq. (|117p depends both on 
orbital phase and frequency, there are truly two coupled, 
first-order differential equations that must be solved. But 
the source term of this equation is incredibly more com- 
plicated than implied here. Even though Sim, Ti„i, Sim 
and pim are known analytically as functions of (j) and 
il, each term contains very long and complicated series 
expansions with fractional exponents that include spe- 
cial functions, such as the polygamma function. For this 
reason, the EOB evolution is not a simple integration, 
as it naively appears to be in this section. Instead, the 
coupled set of first order equations must be solved simul- 
taneously via numerical methods, where here we employ 
a partially optimized, Mathematica routine. 



Disk Modifications to EOB GWs 



1. Effective Hamiltonian 



The effective Hamiltonian controls the conservative 
evolution of the EOB model. We have considered several 
accretion disk effects that directly modify the Hamilto- 
nian, such as the self-gravity of the disk and the increase 
in mass of the SMBH and the CO. Of these effects, the 
latter has been found to be the largest. We concentrate 
on this effect here. 

The increase in the CO's mass can be modeled by solv- 
ing for the time evolution of to*. The differential equa- 
tion that controls this evolution is Eq. where M^^^ 
is given in Eq. (p^ . while ih'^ is given by Eq. (pT|) . No- 
tice that TOj, depends on radius, or on orbital frequency 
by Eq. (|122p . which itself is a function of time. Since the 
accretion rate is not constant, one must solve the system 
of differential equations on Q, and m^, consistently, which 
wc do numerically and perturbatively as follows. First, 
wc solve the frequency evolution equation, setting m^, 
to a constant and neglecting accretion. Second, we use 
Eq. (|I22p to rewrite Eq. in terms of orbital frequency. 
Third, we replace this orbital frequency by the time- 
evolution obtained in vacuum. Fourth, we numerically 
solve the evolution equation for m^,, as its source now 
depends only on time. In step two, we are implicitly dis- 
carding non-linear terms that scale as (m^^Tobs)^, i-C- the 
square of the accretion rate times the observation time. 
This quantity is approximately 10^^ or much smaller 
(clearly much smaller than unity) for typical LISA ob- 
servation times. 

Once the time evolution for the CO's mass has been 
obtained, one must then make sure that all quantities 
that depend on it are properly promoted to time func- 
tions. For example, all quantities that depend on the 
total mass, such as the symmetric mass-ratio rj or the re- 
duced mass /i are modified. In particular, the numerical 
code used to solve the differential equation (Eq. ()II7p ) is 
naturally written in dimensions of the total mass of the 
system, which now becomes a time-function. A simple 
trick to deal with this is to rcscale all mass-quantities by 
the factor M(i)/M(0), which is equal to unity initially, 
but deviates from unity with time. In particular, this im- 
plies that n [M{t)/M{0)]n and f ^ [M{t)/M{0)]-^f. 



GW modeling in the adiabatic EOB framework de- 
pends sensitively on the energy E and the flux J^jw 
Modifications to the MBH or the CO mass naturally 
change all mass scales that depend on the total mass M, 
such as the symmetric mass-ratio 77. Radiation pressure 
and migration modifies the rate of change of the angu- 
lar momentum and thus the fiux. In what follows, we 
explain how wc modify the EOB scheme to account for 
such disk effects. 



Once these substitutions have been made, one can 
solve for the frequency and phase evolutions, with m^,(t) 
a function of time, by providing appropriate initial data. 
When considering BHL accretion, we choose the same 
initial data as in the vacuum case (as explained in the 
end of Sec. lIXBp one year prior to reaching a GW fre- 
quency of 10^^ Hz. The frequency and phase evolution 
can be compared when to* is a constant and when it is 
not, which provides a measure of the effect of BHL ac- 
cretion on GWs. 
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2. Radiation-Reaction Force 

The radiation-reaction force controls the rate at which 
orbits inspiral. This force can be expressed in terms 
of the rate of change of orbital elements, such as bind- 
ing energy, angular momentum and the Carter constant. 
Since we restrict attention to an equatorial, quasi-circular 
EMRI geometry, we need to consider only the energy flux. 

Modifying the EOB model to account for a different 
radiation-reaction force amounts to the rule 

-^GW ^ -^GW ( 1 + ) (123) 

y *Gw J 

in Eq. (jll7p . When modeling an azimuthal wind, 6£ — 
£wmd yjj^ ([SH]). while when modeling migration, then 
S£ is given by Eqs. (|114|) : all throughout £gw is given by 
Eq. (j6|) with e = 0. When substituting in for Si one must 
be careful to use the properly quenched to* (Eq. (ITOp ). 
if the accretion disk effect depends on the rate of BHL 
accretion. All other aspects of the framework can be left 
unchanged, as Eq. (jl23p automatically induces deviations 
from the Kepler relation. 

The system of EOB differential equations can now be 
solved using appropriate initial data. Initial conditions 
for the vacuum and accretion-disk case are prescribed via 
mock evolutions as explained in detail in Sec. IIXBI with 
6£ = and S£ respectively. As explained above, 
this guarantees that both simulations will terminate at 
the same orbital separation. Due to different radiation- 
reaction force laws, however, the starting radii or fre- 
quencies are different in each case for a fixed observation 
time. To account for this, we will later maximize com- 
parison measures over a time and phase shift between 
vacuum and accretion disk waveforms, as we explain in 
SecES 



X. DATA ANALYSIS CONSIDERATIONS 

We are now ready to perform a more detailed data 
analysis study of the accretion disk effects on waveforms. 
We begin by investigating the dephasing of the EOB 
waveforms constructed in the previous section. We then 
continue with an overlap study and end with a discus- 
sion of degeneracies between accretion disk parameters 
and EMRI system parameters. 



A. Dephasing Analysis 

As explained in Sec. Ill B 21 a dephasing study is use- 
ful to roughly determine whether two waveforms can 
be distinguished from each other given a GW detec- 
tion; if the phase difference or dephasing between wave- 
forms is large enough, then they arc distinguishable (see 



Eq. (dni)). Here, we compare the dephasing of the dom- 
inant, (£,to) = (2,2) vacuum and non- vacuum GW 
modes, initialized with the data of Sec. IIX Bl 

If we take one waveform to be "the signal" and the 
other to be "the template," then the dephasing depends 
on two extrinsic parameters contained in the template: 
an overall phase Scf) and time shift St. We here study the 
dephasing after minimizing it with respect to these ex- 
trinsic parameters. The template also depends on other 
parameters, such as the masses and spins, but we here 
set these to be equal to the signal's parameters, i.e. we 
do not minimize the dephasing over such parameters. In 
a realistic data analysis implementation, one would max- 
imize the SNR (or minimize the dephasing) over all pa- 
rameters, at the cost of introducing error into parameter 
estimation due to the use of incorrect templates. 

Before minimizing the dephasing over a time and a 
phase shift, it is worth pointing out that its magnitude 
(computed with relativistic EOB waveforms) roughly 
agrees with the Newtonian results presented in Table H] 
For example, the final dephasing after 12 months of evo- 
lution (initializing the simulations with the same final fre- 
quency and the same initial phase) between vacuum and 
/3-disk migration is ~ 217 rads with the EOB model and 
^ 670 rads with the Newtonian estimates for System I. In 
this case, the relativistic model leads to a larger dephas- 
ing than the Newtonian estimates because, in the latter, 
Sys. I evolves more rapidly due to relativistic corrections 
and the CO spends less time in the Type-II dominated 
region. We have verified that for weakly-relativistic EM- 
RIs, where only Type-I or Type-II migration is in play, 
Newtonian and EOB dephasings agree. 

A more appropriate measure of distinguishability, how- 
ever, requires that one minimizes the dephasing with 
respect to St an d Scj). Following the prescription of 
Eq. (23) in [l65l |. we search for a St and S(j) such that 
\\hit+6t)~f2m < Sf i^nd\\Mt+St)-Mt)-S^\\ < St, 
where 0i^2 and /i^2 sltc the dominant GW phase and fre- 
quencies for waveforms /ii,2(0- The ||j4|| notation stands 
for the integral of A over a time window of length 64Agw , 
where Aqw is the GW wavelength (see e.g. [73 for de- 
tails). The signal hi{t) is assumed to be a vacuum tem- 
plate, while ft,2 is a non- vacuum template for a specific ac- 
cretion disk effect. We choose the tolerances Sf = 10~^^ 
and St — 10~^; decreasing these magnitudes docs not 
visibly change the dephasing results shown below. The 
value of St and Scj) arc unique for a specific set of ft.1.2, 
i.e. for a given accretion disk effect. This alignment pro- 
cedure has been shown to be equivalent to maximizing 
the fitting factor over time and phase of coales cenc e in a 
matched filtering calculation with white noise jl65l | . 

Figure [S] shows the dominant dephasing (left-panel) 
and fractional amplitude difference (right-panel) after 
such alignment. As before, we plot these quantities for 
the dominant GW mode as a function of time in units 
of months, using different color curves for different ac- 
cretion disk effects and different curve styles for different 
types of disks. Observe that after alignment, the dephas- 
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ing increases much less rapidly than in the previous case. 
This implies that the loss of overlap will also grow much 
more slowly. Dephasings after alignment thus correspond 
to the least difference between vacuum and non- vacuum 
waveform phases, without maximizing over intrinsic EOB 
parameters (such as the SMBH's and CO's mass) or other 
extrinsic parameters (such as those associated with the 
observation angles, detector motion, etc.). 

Figure [6] confirms that certain accretion disk effects 
become significant very early during a LISA observation. 
As discussed in Sec. Ill Bl a rough measure of whether 
the dephasing is "distinguishable" for an event with SNR 
~ 10 is whether Stjiaw ^ 1 rad (see Sec. IX Bl for a more 
accurate measure) . The imprint of migration becomes al- 
most immediately distinguishable for Sys. I and /3-disks, 
while it takes at least one full year of observation before 
one can observe the same type of migration for Sys. II 
or BHL accretion for Sys. I. Wind effects also become 
important within one year of observation but for /3-disks 
only. 

Whether an effect is distinguishable is naturally sensi- 
tive to the EMRI parameters and orbital radii. Indeed, 
the bottom panel of Figure [51 corresponding to System- 
II, shows that in this case most effects are much smaller. 
This is mainly due to the assumption that System II's or- 
bit is much closer to the SMBH where disk effects are less 
relevant. The bottom panel further suggests that BHL 
and azimuthal wind effects might be barely distinguish- 
able after 1 year for /3-disks only. Notice that all of these 
findings are consistent with the Newtonian estimates and 
figures presented in the previous sections. 

The right panel of Figure [6] shows the fractional ampli- 
tude difference between the dominant mode of vacuum 
and non-vacuum waveforms after the alignment proce- 
dure described above. The amplitude difference follows 
closely the trend of the dephasing: Type-II migration is 
clearly visible in amplitude changes, while other effects 
are greatly suppressed. The amplitude difference plays 
an important role in the calculation of the overlap and 
the SNR of the difference that we show below. 



B. Overlap Analysis 

Dephasing studies are convenient as rough measure of 
distinguishability for a fixed SNR. However, the SNR 
changes during the GW observation, as signal accumu- 
lates and its frequency enters the detector's more sen- 
sitive domain. A more accurate measure of the detec- 
tion significance of the particular disk effect in the data 
stream is the SNR of the waveform difference between 
the data streams with and without the effect, and the 
so called overlap/mismatch. In Sec. Ill Bl we defined 
all these quantities in terms of cross-correlation integrals 
weighted by the spectral noise of the LISA detector. Note 
that these quantities account for the difference in both 
the phase and amplitude evolution of the GW signal. 

Figure [7] shows the SNR of the difference p{5h) in the 



{£,m) = (2,2) GW harmonic after minimization over 
time and phase shift in units of months. As before, dif- 
ferent curve styles and colors correspond to different ac- 
cretion disk models and effects, as defined in Fig. [TJ The 
thick horizontal line corresponds to an SNR of 10, just 
about the threshold for detection. The numbers at the 
top of this figure show how the vacuum waveform SNR 
builds up as a function of observation time. This figure 
confirms that /3-disk migration is the dominant pertur- 
bation to the measured GW signal and suggests that it 
is distinguishable within 4 months of observation. Mi- 
gration is followed in significance by BHL accretion (for 
either a or /3-disks) and /3-disk azimuthal winds. These 
effects become distinguishable only after a full year of 
observation. All other accretion disk effects arc insignif- 
icant within a 1 year evolution for System I. For System 
II, only /3-disk BHL accretion and azimuthal disks are 
significant and only after 1 year of observation. Other 
effects may become significant for binaries that are closer 
to Earth than 1 Gpc, or if the observation is longer, or if 
the binary orientation relative to the detector is better 
than average. 

This figure is a more realistic estimate of distinguisha- 
bility than the dephasing study presented in the previous 
section. The increase in realism comes at the cost of a 
small drop in distinguishability; e.g. although Fig.[7]sug- 
gests that BHL accretion might be measurable after a 
1 year observation, this is only marginal in the figures 
above. Such a drop is mostly due to the inclusion of 
detector noise in this subsection. Measurability of accre- 
tion disk effect would of course improve if the source is 
closer to Earth, such that the SNR of the signal is larger. 
Irrespective of this, all calculations suggest that Type II 
migration is such a strong effect that it is likely to be 
measurable with LISA. 



C. Degeneracies 

Degeneracies between EMRI system parameters, such 
as the SMBH and CO's mass, and accretion disk param- 
eters could deteriorate the extraction of accretion disk 
parameters from EMRI GW observations. If one were to 
maximize the overlap function over all parameters (in- 
stead of just a time and phase offset, as done in the pre- 
vious section), one might find mismatches much closer 
to zero, at the cost of biasing parameter extraction. In 
this subsection, we investigate this issue and the spectral 
signature left by disk-induced effects. 

The effect of possible degeneracies can be assessed by 
investigating the Fourier transform of the GW response 
function, as this is the main ingredient in matched fil- 
tering. We restrict our study of degeneracies to a simple 
analytical estimate of the Fourier transform using the 
New tonian st ationary phase approximation (SPA) (see 
e.g. |l70l . llTll ). First, let us review this approximation 
in vacuum GR, and then consider the modifications in- 
troduced by leading-order disk effects. 
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FIG. 6. Aligned dephasing (left) and aligned fractional amplitude difference (right) as a function of time in units of months 
for the dominant GW mode. Line style and color follow the same notation as in Fig. [5] The thick solid lines signal a 1 radian 
dephasing. The top and bottom panels correspond to System I and II respectively. Observe that the minimized dephasing 
exceeds unity in a short observation time for /?-disk migration and System I, while for System II, only BHL accretion and 
/3-disk wind effects do so after a full-year of observation. 
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FIG. 7. SNR of the aligned, dominant waveform difference as a function of time in units of months. The left panel corresponds 
to System I and the right one to System II. Different curve colors and styles correspond to different disk effects and disk models, 
respectively, as in Fig. [5] The thick horizontal dashed line correspond to an SNR of 10. Observe that the SNR of the difference 
exceeds the SNR threshold in four months for /3 disk migration and Sys. I, while it take a full year of observation for the same 
to occur when modeling BHL accretion and /?-disk wind effects. 



Tlic Fourier transform of tlie response function h{t) 



ml) 



h{f)= / h{t)e'^'^'dt. 



(124) 




i{2irfto-4>o) 



(125) 



Tills generalized Fourier integral can be solved via tlie 
method of steepest descent, assuming the amplitude 
changes slowly relative to the phase and noting that the 
complex phase ip = iirft — (/)gw has a stationary point at 
d'ip{f,to)/dt = 2Trf — d(f>Q^(tQ) / dt = 0. In this approxi- 
mation, the Fourier transform becomes (see e.g. Eq. (4.5) 



where the factor of 8/5 accounts for sky-averaging 
over beam pattern functions. The quantities 
[tQ,\p{f,to),il){f,to)\ for a fixed / can be found by 
assuming that the phase and time of merger are fixed 
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where r^o = //2 is the stationary point and fo = f(//2). 
The quantity f can be constructed from {dE/df)~^E 
and E = J^awir,^) is the GW energy loss rate, 
given in Eqs. (jlisp in the Newtonian approximation or 
Eqs. (|118lll9p in the EOB approximation, if we neglect 
accretion disk effects. 

Neglecting disk effects, we can readily evaluate each 
of these terms to leading (Newtonian) order. Using that 
the orbital frequency is given by Eq. 1^ and the rate of 
change of the orbital separation by Eq. Eq. (|126p 
becomes to ^ tc - i5/256)Mu-^/^, while Eq. 07\ be- 
comes (j){to) = (l)c- (l/16)u-5/3, where we have defined 
the reduced frequency u = nMf, with the chirp mass 
Ai = g^/^M.and q = m^^/AI, the mass-ratio. Using the 
Newtonian expr essio ns for A{t) — M/Dl{M^Y/^ (see 
e.g. Eq.(3.5) in [ml), the sky -averaged Fourier ampli- 
tude becomes 
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where is the luminosity distance of the source and / 
is the observed GW frequency, while the Fourier phase is 



^^vac(to,/) = 27r/to-0(/) 



-5/3 



128 



const, (129) 



Let us now repeat this calculation with accretion disk 
modifications. We modify the above algorithm by re- 
placing E = i?vac(l + S£/£gw), where 6£ is given by 
Eqs. (Unni), (UnS) or dnni) for Type I, H-SC and II-IPP 
migration respectively, or Eq. (|80p for azimuthal winds. 
We replaced to* by Eq. ([TTjl when modeling unquenched, 
BHL accretion. The resulting frequency-domain phase 
and amplitude can be parameterized as 



and 



l^l/l^lvac = 1 - Bial^m, 



C2 ]\,f°-3„<i4 



(130) 



(131) 



where qo = q/10^^ is the normalized mass-ratio q — 
M^/nii, and uq = (7rA^/)/(6.15 x 10~^) is a normalized 
reduced frequency and a GW frequency of 10"^ Hz. The 
parameters {Ai, Bi^hi) are given in Table Hill while no- 
tice that (ci, C2) are the same as those in Table HH Note 
that the expressions in Eqs. (|130p and (|13ip are valid 
only in the regime of frequency space where the accre- 
tion disk effects are small perturbations away from the 
vacuum evolution (ie. at sufficiently small separations). 
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TABLE III. Columns are parameters in Eq. (|130|) and rows are 
migration effects. As in Table U the notation x (y) = x x 10^ 
in radians for Ai and dimensionless for Bi . Observe that the 
frequency exponent 0.5 < 0, implying that these accretion disk 
effects are dominant at small frequencies (large radii). 



Let us discuss these results further. First, notice that 
corrections to "0(^0, /) due to Type II migration are or- 
ders of magnitude larger than all other effects, as shown 
by the magnitude of Ai. Second, notice that all disk- 
induced corrections depend on negative powers of fre- 
quency (or reduced frequency u in this case). This is be- 
cause such accretion disk corrections are largest for large 
radii, equivalent to weak-field GR effects. In fact, they 
are dominant over the leading-order vacuum term (the 
factor of u~^/3 i/vac) at low frequency. This suggests 
that migration effects are not strongly correlated with 
GR vacuum terms in the PN approximatiorP^I. as the 
latter depend on positive powers of u relative to u~^/^. 

One might wonder how the accretion disk effects mod- 
ify the Fourier phase and amplitude when they are 
not necessary a small perturbation away from the vac- 
uum evolution. In general, the accretion disk correction 
changes the functional form of the phase or amplitude as 
follows: 



l + Aar<fAC|g° 



(132) 



where j/vac = ("0, 1'^l) when A = {Ai,Bi). This means 
that, unlike what Eqs. (|130p and (|13ip suggest, the ac- 
cretion disk effects always suppress the vacuum evolution 
as A > 0. Figure [5] shows the absolute value squared of 
the Fourier amplitudes as a function of frequency for an 
EMRI with M, = IO^Mq and to* = IOMq and dif- 
ferent accretion disk effects (neglecting all quenching). 
For comparison, we also plot the vacuum amplitude and 
the spectral noise density curve. Observe that below 
/ < 10"'^ Hz, accretion disk induced migration becomes 
dominant over GW emission, and the Fourier amplitude 



However, the detected signal is also modulated during LISA'S 
orbit around the Sun which we neglect in this paper. This can 
introduce correlations between parameters that change the GW 
phase slowly, su ch as migration effects, source direction, and ori- 
entation angles Il72ll . 



34 





a Gap 


1 M 1 1 1 II 

PGap 


ISCO 

/ 


1 # 

/ 
















: 












vacuum 


■ 




X f ' 








BHL a 






/ 1 ■ ''/X .'-'^^ 








BHL p 


- 




/ ^ I •!* 








Wa 






/ / \ ' ' \ 
I • f /• .'1 








wp 






/ ■ '7 "1 








Ml a 






/ > / / ■" I 








Ml p 






// .* / .■ I 








M2aj.p 




-/a 
























M2P,^ 














M2P,„ 






















M.I < < . 




III 







10"^ 10"** 10"'' 10"" 



f[Hz] 

FIG. 8. Absolute value squared of the Fourier amplitudes 
with different accretion disk effects (see label) as a function 
of frequency in units of Hertz. We also plot this amplitude in 
vacuum and LISA's spectral noise density. Vertical lines from 
right to left correspond to the ISCO frequency, the frequency 
at which a gap forms for /3-disks and for a-disks. 



is significantly different. This effect was also demon- 
strated todecrease the GW background for pulsar timing 
arrays |l48l | . A gap is expected to be opened at small fre- 
quencies, but to close at the radii indicated with vertical 
lines (see Sec lIIIBj) . At these frequencies Type-II migra- 
tion transitions to Type-I. Coincidentally, hydrodynamic 
drag due to an azimuthal wind and BHL accretion cease 
at smaller frequencies. 

The precise accuracies to which disk parameters can be 
estimated is difficult to ascertain. A crude Fisher analy- 
sis that neglects degeneracies (a diagonal approximation) 
suggests extraction accuracies of up to 1%. We expect, 
however, that correlations will deteriorate the accuracy 
of extraction down to 10% [t^. Ultimately, a proper as- 
sessment of the accuracy to which disk parameters could 
be extracted from EMRI observations would require the 
detailed mapping of the likelihood surface with relativis- 
tic EMRI signals, full Fourier transforms, and improved 
disk modeling (including relativistic effects and magnetic 
fields). This, however, is beyond the scope of this paper. 



XI. DISCUSSION 

We have explored the effect of accretion disks on the 
GWs emitted during the inspiral of a small BH into a 
much more massive one. We have found that disk migra- 
tion has the biggest influence on the EMRI dynamics. If 
EMRIs are detected with LISA, our study suggests that 
migration could be measurable within 4 months of obser- 
vation for a /3-disk. Depending on the particular EMRI 
considered, a gap could open in the accretion disk, leav- 
ing an imprint in the GW observable. We have studied 
possible degeneracies between accretion disk and EMRI 
parameters and found that they are weakly correlated. 



This is because disk effects are dominant at large sep- 
arations, thus introducing "negative" PN terms in the 
frequency-domain waveform phase. The LISA detection 
of a GW from an EMRI embedded in an accretion disk 
could therefore allow us to measure a combination of the 
a and rii, parameters of the disk. 

The prospect of probing accretion disks with GWs can 
be improved in the presence of an EM counterpart. If the 
EM luminosity of the SMBH's accretion disk is observed, 
one might be able to separately measure both a and m, 
by adding the information from a GW detection. If the 
GW detection does not show evidence of accretion disk 
effects but an EM signal is present, then one might be 
able to distinguish between a and /3-disks. 

Several caveats should be kept in mind about elements 
of our analysis that could be improved on in the future. 
First, we have considered a very specific type of EMRI, 
consisting of quasi-circular orbits on the equatorial plane. 
In principle, the CO could be in an inclined and eccen- 
tric orbit. If so, certain accretion disk effects that were 
negligible in our analysis, such as the self-gravity of the 
disk, could become more important. Such effects would 
induce additional apsidal and nodal precession, which 
could leave a detectable signature on both the EM and 
GW signals. 

Another important caveat is that we did not explore 
possible degeneracies with other EMRI parameters, like 
source direction. It is possible that the orbital modula- 
tion of LISA can induces a time dependence mimicing 
the effects of accretion disks. However, since these mod- 
ulations are periodic, we expect no strong degeneracies 
after a multi-year observation. 

Another important issue involves the accuracy of the 
EMRI model used. Although recently an EOB-inspired 
EMRI waveform model was developed [tsI - ItB} . no such 
model has been thoroughly tested for generic EMRIs. 
Even in the case of quasi-circular inspirals on the equa- 
torial plane, the FOB model is still only accurate up to 
dephasings of order 10 radians. This implies that ac- 
cretion disk effects associated with dephasing signatures 
of this magnitude, such as the effect of BHL accretion, 
might be difficult to disentangle. However, accretion disk 
effects are also important at large separations, whereas 
mismodeling errors are of concern only close to the ISCO. 
Thus, the detection of EMRIs at moderate binary sepa- 
rations might allow the extraction of small accretion disk 
effects. Moreover, the effects of migration are much larger 
than any possible waveform mismodeling and so they can 
be easily isolated. Ultimately, a more detailed Markov- 
Chain Monte-Carlo study is required to determine the 
accuracy to which disk effects can be measured. 

Finally, let us highlight the large uncertainty that cur- 
rently exists in accretion disk modeling. Even when con- 
sidering radiatively efficient thin Newtonian disks, there 
are several viable models, here parametrized as a and /3- 
disks, that can have drastically different effects on EMRI 
inspirals. Unfortunately, these Newtonian models are ex- 
pected to be highly inaccurate precisely in the regime 
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where EMRIs are most easily seen, i.e. in the relativistic 
regime close to the SMBH. A natural extension of this 
work would be to consid er EMRI s in the context of rel- 
ativistic accretion disks |lOll Il02| , for exaniple through 
the incl usion of r elativistic BHL accretion [381 and mi- 
gration plITiii. 



ACKNOWLEDGMENTS 

We are grateful to Cole Miller, Ramesh Narayan, 
Shane Davis, Scott Hughes, and Zoltan Haiman for use- 
ful comments and suggestions. BK and NY acknowl- 
edge support from the NASA through Einstein Postdoc- 
toral Fellowship Award Number PF9-00063 and PFO- 
110080 issued by the Chandra X-ray Observatory Center, 
which is operated by the Smithsonian Astrophysical Ob- 
servatory for and on behalf of the National Aeronautics 
Space Administration under contract NAS8-03060. BK 
acknowledges partial support by OTKA grant 68228. AL 
acknowledges support from NSF grant AST-0907890 and 
NASA grants NNX08AL43G and NNA09DB30A. 



Appendix A: Phase Correction due to CO Accretion 



where 



ro fro 
df 



(A5) 



Substituting in Eq. (|A2[) , the initial separation becomes 




i^m*)3,3 



4096 mi ^T 



1/4 



(A6) 

where t^pa is the dimensionless observation time given by 
Eq. the m index was introduced to distinguish from 
other modifications below, c.f. Eq. for 7~o. Thus, the 
change in the lower integration bound in Eq. (jAip is 



5f™* = - fo 



5 AP 



Since the total relative change in is very small 
during the observation, we may approximate 1/771^,(1) in 
Eq. (|XT|) as 



1 



1 



so that Eq. (jAip becomes 



2 

*,0 "'*( 



m^,{t')dt' 
(A8) 



Let us compute here the GW phase evolution when 
the CO's mass is varying due to BHL accretion. The 
GW phase is given by the first line of Eq. ([H)) . 

ff rif 5 ri f3/2 

r,,, r 32 Jf^^, m^r) 

Two important corrections arc induced by the CO's vari- 
able mass: a change in the limits of integration fo — >■ fo' 
and a change in the denominator of the integrand in the 
phase evolution, to* — ^ m*(t). We expand both in a 
Taylor series around fg ^-nd the initial CO mass to*,o at 
f = fo' . We evaluate the change in the total GW phase 
induced by accretion while keeping ff, the observation 
time T ~ tf ~ Iq, and m^^^o fixed. 

The limits of integration can be computed by integrat- 
ing the radial inspiral evolution equation ([8]) with the 
time dependent m^(t), 

,-4_4^_^_L£^^(t)dt. (A2) 

where 
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V- _ f3/2 
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1 M. 



16 TO 



/2 ^5/2 



25 M, 



2048 TO^ o 



r~ ('^"^*)3/2,3 ' 

(A9) 



where the to,^ label denotes that the CO is accreting, 
in the second line we have changed integration variables 
using df = df' jf' and used Eq. (jA5p . 

Relative to the GW phase without accretion, 



32to^,o ° 
25 Ml ( {Sm^) 



25 M3 

— (OTO*) 



2048 to2 



3/2,3 



3,3 



2048 



-3/2 



'3/2,3 



(AlO) 



Equation (jAlOp and (jA5P can be used the calculate the 
GW phase shift for arbitrary to^,. In the next two sub- 
sections we consider BHL accretion and gas supply 
limited BHL accretion = M,. 



Jt„, 



(A3) 



Note that rh^ = TO*(f) given by Eqs. (|54l70p . Changing 
integration variable from t to f using dt = df/f using 
Eq. ^ for r(f), we get 



1. Unsaturated BHL accretion 

First consider the case where BHL accretion is not lim- 
ited by the amount of local gas supply, m^{r) < M,. 



25 Af^ \ 
(<5to.)3_3 I , (A4) 



4096 m 



*.o 



Since ((5m<.)^ g <JC mi,fiT, we approximate the lower integration 
bounds in Eq. I IASI I with fg. 
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BHL accretion is not quenched by local gas supply if the 
observation is limited to orbital radii r < fq, or equiva- 
lently, if the dimensionless accretion rate onto the SMBH 
TO, > TO,_q. In practice, this is the case during a year of 
observation approaching ISCO, if the SMBH accretion 
rate is moderate to high rn, > 0.3 or if the mass-ratio is 
very small to*/M, < 10~^ [see Eq. ((M)) ]. 

Substituting the BHL accretion rate, Eqs. (|54IA5p for 
a and /3-disks into Eq. (|A10p yields 
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(A12) 



where to*i refers to the initial CO mass at ro, x = r//ro, 
^0 = roirf,T) is given by Eq. ([T I|, and fp 20 = ro/20. 
The phase evolution given by Eq. (jAlH - ET2]) depends on 
only two sets of parameters: the time-independent coeffi- 
cient preceding the parentheses, and the time dependent 
quantity T/fj. The later also appears in the standard 
inspiral phase expression, Eq. (jl4p . which we used to dis- 
tinguish two cases: when the observation time is long or 
short relative to the inspiral timescale at the given radius 
[see Eqs. ([15]-[l6l)]. 

We can similarly distinguish here between two asymp- 
totic cases. For long observations (T ^ Tcrit) or small 



separations (r/ ^ ^/,crit), ''o ~ ^ 



1/4 



-1/4 



« 1, 



we can approximate the dephasing with the general for- 
mula 



^GW.long 
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rpCe 

± yr 



(A13) 

where the coefficients (Ci,Cj) are given by the first two 
rows of Table |TT] and n = 2 for Type H migration, while 
n = 1 for all other migration disk effects. For short ob- 
servations (T ^ Tcrit) or small separations (f/ ^ ^/^crit), 
r/r^ <C 1, fp f/ -|- r/(4f^), a; « 1 — T/(4fj), and we 
can approximate the dephasing with the general formula 



^GW, short 



I yr'/io 



1 -£). 



Ml 



5 '/,10/ 

(A14) 

where the coefficients {Di,dj) are given by the first two 
rows of Table HI] and n = 2 for Type H migration, while 
n = 1 for all other migration disk effects. 



2. Gas supply limited BHL accretion 

If TO^(f) > M,, the accretion onto the CO is limited by 
the amount of local gas supply. This is the case outside 
f > fq, or if the SMBH accretion rate satisfies to, > m,,q, 
see Eq. for (fq,TO,_q). In practice, this is the case 
for TO, < 0.1 for beta disks if the final orbital radius 
is not very close to the ISCO (e.g. f/ > 10). For larger 
TO, accretion rates approaching the ISCO or intermediate 
mass-ratio inspirals, BHL accretion starts supply limited 
and becomes unsaturated near the ISCO. 

Assuming to* = A/, in Eqs. (jASIAlOp . we get 
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'0,20 
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(A15) 



for both a and /3-disks, where again fp = fQ{ff,T) given 
by Eq. ([TT|) and x = f//fp. Relative to the BHL accre- 
tion case with unlimited gas supply Eq. (jAlllET2)) . the 
phase is a less steep function of radius. The evolution is 
again determined by two combination of parameters, the 
constant coefficient (which is now different) and the time 
dependent quantity T/fj. 

In the two limiting cases, where the observation is long 
or short relative to the inspiral timescale, we again can 
parameterize the dephasing as in Eq. (|AT3l) and (|AT4l) 
respectively, where the coefficients {Ci,Cj) and {Di,dj) 
are given by the third row of Table HIl 



Appendix B: Phase Correction due to Wind and 
Migration 

In this appendix wc derive the correction to the GW 
phase due to a modification in angular momentum dissi- 
pation. The result of a similar, but more general, angular 
momentum dissipation rate = -^f^ is presented 

in Paper I [s^. Wc focus on quasi-circular orbits only. 
From Eq. (O, 

2 2 ro(.)(^?)'7s)* 



dt 



dr 



'■^ , df 



(Bl) 



where we have expressed the GW phase as a function 
of radial and temporal derivatives of the specific angular 
momentum of the CO. Similar to Appendix we calcu- 
late the total total change in the GW phase relative to 
the unperturbed GW inspiral phase, by keeping the final 
separation f j and the observation time T fixed. Relative 
to the vacuum inspiral Eq. (|14p . a modified angular mo- 
mentum dissipation causes a phase shift by changing the 
lower integration bound and I in Eq. (|Bip . Since the ad- 
ditional specific angular momentum loss hi is small rela- 
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tive to ^Gw, we may approximate the result by expanding 
in a series in the small quantity 8lj 1^^. 

First to estimate fp' , note that the accelerated dissipa- 
tion of angular momentum causes an accelerated inspiral 
rate. For circular orbits, I = M,f^^^, E = Af*/(2f), so 
that 



2fi/2 ^ 
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M. 5 M2 \^ I 

This equation can be integrated to give 
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(B2) 



(B3) 



where we have introduce the index i to distinguish from 
other effects, and wc dcfinco 



dr -. — dr' . 



(B4) 



Comparing to Eq. the change in the lower integra- 
tion bound 



(B5) 



Now using 



1. Wind — unsaturated BHL accretion 

If the additional angular momentum dissipation is 
caused by an azimuthal wind where the accretion onto 
the CO is the unsaturated BHL rate, 5i{f) = (^^i^^^{r) 



given by Eq. 
we get 



81]). Substituting in Eqs. dBg]) and ([14 
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where x = rf/fo and ro = ^o(ffTT), see discussion fol- 
lowing Eqs. (|A12IA15[) . Following the derivation pre- 
sented in Appendix I A 11 wc can expand in long and short 
observation times relative to the inspiral timescale to re- 
cover a dcphasing as in Eq. (KT3\i and (|JT4)) with the 
coefficients (Ci^Cj) and (Di,dj) given by the fourth and 
fifth rows of Table m 
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for the integrand in Eq. (jBl)) . the GW phase is 
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so that relative to the vacuum inspiral phase 
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Equations ((B8)) and (|B4| are the analogues of Eqs (|A10|) 
and (jASP for the GW phase shift caused by an addi- 
tional source of angular momentum dissipation Si. In the 
next three subsections we consider the azimuthal wind, 
Type-I migration, and Type-II migration. In these cases 
Si/iaw = -^f^ where A and B are constants, so the 
integral in Eq. (|B4p can be evaluated analytically. 



2. Wind — quenched BHL accretion 

If the additional angular momentum dissipation is 
caused by an azimuthal wind where the accretion onto 
the CO is limited by the amount of gas supply then we 
substitute Siif ) = Ci^-i^"^(^) Eq. (EH) mto Eq. ^ to 
obtain the corresponding total GW phase shift 
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where throughout this section the top and bottom rows 
correspond to a and /3-disks, respectively, and the pa- 
rameters X = x{ff,T) and 7~o = fo(r/,T), see discussion 
following Eqs. (jA12IA15p . Notice that although we are 
here setting m^, = M, , there is still a dependance on the 
type of accretion disk, due to the factor of 7 in Eq. ([50]) . 

Again, as in Appendix IA 11 we expand in long and short 
observation times to obtain the dephasing of Eqs. (|A13p 
and (|A14p with the coefficients (Ci, Cj) and {Di, dj) given 
by the sixth and seventh rows of Table [TTl 



3. Type-I migration 



Since ^wind/^Gw <SC 1, we approximate the lower integral: 
bound with fg in Eq. I IB4II . 



iion Let us now compute the dephasing introduced by the 
angular momentum dissipation due to Typc-I migration. 
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Substituting Eqs. (|106|) into Eqs. (|B8|) and (|B4p we get given by 
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,/3,TI _ _i -4/5 "^.1 42/5 
"gw — -l^-O "1 '0,20 
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where again x = x{ff,T) and fo = ^0(^/7 

We now take the two Hmiting cases of long and short 
observations. For long and short observations, we find 
dephasings as in Eq. (jA13p and (|A14p with coefficients 
{Ci,Cj) and {Di,dj) given by the eighth and ninth rows 
of Table El 



4. Type-II migration 

Consider now the dephasing corresponding to the an- 
gular momentum dissipation due to Type-II migration. 
Substituting Eqs. (fT08)) and ([TTO| into Eqs. and 
(|B4| we get 
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where now both the angular velocity fi(f ) and the inspiral 
rate f are modified by the disk potential, sec Eq. (P7| and 
(|102p . The later also implies a change in the integration 
bound for a fixed ff and observation time. 

We can integrate Eq. (|102p perturbatively to obtain 
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where we have defined 
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Substituting Eqs. ^ and in Eq. JCT]) gives 
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Subtracting from this expression the vacuum expression 
for the GW phase, we find 
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where again x = x{ff,T) and fo = fo(f/,T). No- 
tice that we have normalized the CO's mass to m^,■2 — 
m^,/(100Af,), where a gap opens and Type-II migration 
occurs. 

Let us now take the two limiting cases of long and 
short observations. For long and short observations, we 
find a dephasing as in Eq. (|A13p and (jA14p respectively, 
with {Ci,Cj) and {Di,dj) given by the tenth, eleventh 
and twelfth rows of Table |lTl 



Appendix C: Phase Correction due to the Disk 
Gravity 

Let us study how the gravitational potential generated 
by the disk, $disk, affects the GW phase. The latter is 



idisk 
''gw 
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(G5) 

The result is algebraically similar to the phase shift due 
to a modified angular momentum loss rate Eq. (jB8p or 
a modified CO mass Eq. (jAlOp . Note that the overall 
minus sign is compensated for by the sign of the potential 

«'disk < 0. 



1. Disk Potential 

Let us compute a more convenient form of the disk 
potential $disk- Assuming that = Sq?"''', the accre- 
tion disk parameters (a,TO,,M,), carried by Eg, can be 
taken out of the integrals in Eqs. (|91l92p . The integral 
only depends on (7, f , fmin, ?^max)- Wc may evaluate the 
integrals in the Legendre expansion (|92p . 
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Next we exercise the gauge freedom to set 
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and in the following drop the "new" specifier. Note that 
for both disk models 7 > —1; it is 3/2 and —3/5 for a. and 
/3-disks. The disk potential in Eq. (|C6|) is expressed as a 
sum of two terms in the eurly brackets, which correspond 
to the potential of the disk interior and exterior to the 
orbit, respectively. 

Let us discuss the contribution of various multipolar 
harmonics, i. 
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In each case we have arranged the terms in increasing or- 
der for 7 = 3/2. Equation (|C8[) shows that for 7 > 1 and 
r <C fmax the potential is dominated by the quadrupolar 
harmonic 1 — 1^ 



"2(7-1) • 
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where 
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1 r(i + i)r(V) 



7- 



1 



r(i^)r( 



(Cll) 



so that Co = 1.38 for /3-disks where 7 = —3/5. Equa- 
tions (jC9[) and (jClOP represent the asymptotic solutions 
for a and /3-disks, respectively. Substituting the particu- 
lar surface density profiles for the two disk models leads 
to the potential given by Eq. (jHS)). 



2. Axisymmetric disk gravity without a gap 

We can now evaluate the GW phase shift, Eq. (jC5p . for 
the potential generated by a and /3-disks without a gap. 
We restrict to fmin '^f <^ fmax and substitute Eq. ([95]) 



5 ^ -20/21 ^^.5 -11/2 
' n 



13/21 '0,20 



10 ^ 10 J ' 



6^'j8'^ = 3.0xl0-^ai 



. 3/5 , .11/5 

-4 -4/5 ^,2 iVj,g _39/10 
"0,20 



13 



10 



^39/10 _ if y27/5 

10' 



(C12) 



(C13) 



When — 1 < 7 < 1, the disk potential is asymptotically 
independent of the inner and outer boundaries for f,„i,i <C 
^ <SC ^max ■ We can then analytically evaluate the infinite 
sum, and get 



$^<i(f)«27rcoSof^+i 



(CIO) 



where again x — ff/vQ. The sign difference is due to the 
fact that the disk exerts an outward pull for a and an 
inward push for /3-disks. In the long and short observa- 
tion limits, we find dephasing as in Eq. (jA13[) and (|A14[) 
with {Ci,Cj) and {Di,dj) given by the thirteenth and 
fourteenth rows of Table |TT1 
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